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I.  STATEMENT  OF  THE  PROBLEM 


Consider  the  microstrip  discontinuity  shown  in  Figs.  1  and  2.  In 
region  a  (-00  <  z  <  z )  of  Fig.  1,  an  infinitesimally  thin  perfectly  con¬ 

ducting  strip  of  width  w  and  height  h  above  a  perfectly  conducting 

3  3 

ground  plane  runs  parallel  to  the  z  axis.  In  region  c  (zc  <  z  <  <*>)  of 
Fig.  1,  an  infinitesimally  thin  perfectly  conducting  strip  of  width  w 

c 

and  height  h^  above  the  ground  plane  runs  parallel  to  the  z  axis.  These 

two  strips  are  called  strips  a  and  c.  They  communicate  by  means  of  what 

is  called  a  via  in  region  b  ( z £  z  _<  Z£)  of  Fig.  1.  It  is  assumed  that 

the  extent  (z  -z  )  of  the  via  in  the  z  direction  is  greater  than  or  equal 
c  a 

to  zero.  The  via  is  a  perfectly  conducting  surface  of  arbitrary  shape. 

It  is  assumed  that  the  via  establishes  an  electrical  connection  between 
the  two  strips. 

Above  the  ground  plane,  the  strips  and  via  are  immersed  in  a 
medium  that  has  constant  permeability  y.  Although  both  of  the  strips 
run  in  the  z  direction,  they  may  be  offset  from  each  other  in  the  trans¬ 
verse  plane  as  shown  in  Fig.  2.  In  Fig.  2,  the  strip  in  region  a  begins 
at  (x**a)  while  that  in  region  c  begins  at  (x*xc>. 

The  objective  is  to  calculate  the  excess  inductance  of  the  struc¬ 
ture  shown  in  Fig.  1.  This  excess  inductance  is  called  Lg  and  is  de¬ 
fined  by 


L  -  lim 

e  £  -*C0 

a 

£  -*» 
c 


\p  -  £  ili 

/  a  a 


£  ip 
c  c, 


(1) 


In  (1),  the  quantities  ip,  ipg,  and  \pc  are  defined  in  terms  of  the 


I 


Fig.  1.  Two  perfectly  conducting  strips  above  a  ground  plane. 

strip  in  region  a  communicates  with  that  in  region  c 
means  of  the  via  in  region  b. 
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Cross  sectional  view  of  the  two  strips  and  the  ground 
plane.  The  via  is  not  shown  here. 


magnetostatic  problem  in  which  flu*  total  current  I  flows  i rom  (z  =  -  -a) 

on  strip  a  through  the  via  to  (z  3  '")  on  strip  r. 

The  quantity  if-  is  the  magnetic  flux  which  passes  through  a  loop 

formed  by  the  portion  of  the  strips  and  via  for  which  (-2  ‘  z  <  2 ■  ) , 

a  —  —  c. 

the  ground  plane,  and  connecting  lines  at  (z  =  -  V.  )  and  at  (z  =  ii^). 

The  reference  direction  for  ig  is  into  the  paper  in  Fig.  1. 

The  quantity  t p  is  the  magnetic  flux  per  unit  length  of  strip  a 

cl 

far  from  the  via.  More  precisely,  if  z  is  such  that  the  intersection  of 

o 

strip  a  with  the  (z=z  )  plane  is  far  from  the  via,  then  ij>  is  the  magnetic 

flux  which  passes  through  a  loop  formed  by  the  portion  of  strip  a  for  which 

(zq  -  ’  m  _  z  •'  ZQ),  the  ground  plane,  and  connecting  lines  at  (z  ■  zq  -  lm) 

and  at  (z  =  z  ) .  Here,  "lm"  denotes  one  meter  and  represents  a  unit  length 
o 

along  the  z  axis.  The  reference  direction  for  ip_  is  into  the  paper  in 

cl 

Fig.  1. 

The  quantity  tp^  is  the  magnetic  flux  per  unit  length  of  strip  c  far 

from  the  via.  More  precisely,  if  z  is  such  that  the  intersection  of  strip  c 

o 

with  the  (z-z^)  plane  is  far  from  the  via,  then  <P^  is  the  magnetic  flux  which 

passes  through  a  loop  formed  by  the  portion  of  strip  c  for  which 

(z  •'  z  '  z  +  lm),  the  ground  plane,  and  connecting  lines  at  (z=z  )  and  at 
o  —  —  o  ^ 

(z=z  +lm) .  The  reference  direction  for  i|)  is  into  the  paper  in  Fig.  1. 
o  c 

Because  alL  three  quantities  <p,  ip  ,  and  i|/  are  proportional  to  1, 

3  C 

expression  (1)  for  L  does  not  depend  on  1 .  If  v  is  sufficiently  large, 

e  a 

the  addition  of  any  positive  quantity  A*  to  1  increases  \p  by  A 9  so 

that  the  quantity  (p  -  f.  ip  )  on  the  right-hand  side  of  (1)  does  not  change. 

3  3 

If  V.  is  sufficiently  large,  the  addition  of  any  positive  quantity  tJ.  to 
c  c 

9.  increases  p  by  !-.9.  i|>  so  that  the  quantity  (ip  -  i  j  )  on  the  right-hand 

M  '  J  f*  0  C  C 


side  of  (1)  does  not  change.  Therefore,  the  limit  In  (1)  exists  and  is 
Independent  of  I.  Expression  (1)  for  Lg  is  similar  to  [1,  Eq.  (1)]  for 
the  excess  inductance  of  a  microstrip  right-angle  bend. 


If  the  origin  of  the  z  coordinate  were  shifted  a  distance  A*  to  the 

right  and  if  the  loop  through  which  the  flux  threads  were  physically  the 

same,  then  would  not  change,  H  would  increase  by  Az,  and  l  would  de- 

a  c 

crease  by  Az.  Therefore,  the  right-hand  side  of  (1)  would  increase  by  the 

amount  Az(i|>  -  4>  ) /I.  Thus,  as  defined  by  (1),  L  depends  on  the  origin  of 
x  "a  c 

the  z  coordinate. 


II.  EXCESS  INDUCTANCE  IN  TERMS  OF  MAGNETOSTATIC  CURRENTS 

In  this  section,  the  limit  on  the  right-hand  side  of  (1)  is  ex¬ 
pressed  in  terms  of  magnetostatic  currents.  Consider  the  right-hand  aide 
of  (1).  Since  all  three  quantities  ip*  4  jl>  »  and  approach  infinity  as 

SL  and  l  approach  infinity,  the  two  subtractions  in  (1)  must  be  performed 
&  c 

before  passing  to  the  limit.  The  flux  $  is  given  by 


* 


l  •  r 


ds 


(2) 


where  is  the  magnetostatic  current  density  on  the  entirety  (-00  <  z  <  °°) 
of  the  strips  and  via.  The  amplitude  of  J  is  adjusted  so  that  the  total 
current  flowing  from  left  to  right  on  the  strips  and  via  in  Fig.  1  is  I. 
In  (2),  j$(J)  is  the  magnetic  field  due  to  the  existence  of  .J  in  the  pre¬ 
sence  of  the  ground  plane,  is  a  surface  that  caps  the  loop  described 
in  the  paragraph  that  follows  (1),  n  is  a  unit  vector  normal  to  S^,  and 


6 


ds  is  the  differential  element  of  surface  area.  The  vector  n  points 
into  the  paper  in  Fig.  1. 

In  the  following  manner,  the  integral  over  in  (2)  is  trans¬ 
formed  into  an  integral  over  C  where  C  is  a  contour  that  goes  from 
(z  *  -fca)  on  strip  a  through  the  via  to  (z“^c)  on  strip  c.  First,  B(J) 
is  expressed  as 


B(J)  =  V  x  A(J) 


(3) 


where  A(J)  is  the  magnetic  vector  potential  due  to  the  existence  of  J 
in  the  presence  of  the  ground  plane.  Substitution  of  (3)  into  (2)  yields 


*  = 


(V  x  A(J))  •  n  ds 


(4) 


Stokes'  theorem  [2,  Eq.  (42)  on  p.  489]  reduces  (4)  to 


* 


J  A(J)  •  dfc 


(5) 


where  the  contour  is  the  loop  described  in  the  paragraph  that 
follows  (1).  On  the  strips  and  via,  proceeds  from  left  to  right  in 

Fig.  1.  In  (5),  d®.  is  the  differential  of  the  radius  vector  on  C^.  From 
the  method  of  images,  we  have 


A  tan  (J)  =  0  011  the  ground  plane  (6) 

where  the  subscript  tan  denotes  the  component  tangent  to  the  ground 
plane.  According  to  (6),  the  part  of  C{J  on  the  ground  plane  does  not 
contribute  to  the  integral  in  (5).  Far  from  the  via,  J  is  entirely  z 
directed  and  so  is  A.  Therefore,  if  ®  and  9  are  sufficiently  large, 

3  C 


the  parts  of  C„  along  the  connecting  lines  at  (z  -  -l  )  and  at  (z  ■  &  ) 
do  not  contribute  to  the  integral  in  (5).  As  a  result,  (5)  reduces  to 

' 

ip  «  A(J)  •  d£  (7) 

C 

where  C  is  defined  in  the  first  sentence  of  this  paragraph. 

The  magnetostatic  current  density  J  exists  on  S  where  S  is 

-  CO  CO 

the  entire  surface  (~°°  <  z  <  00 )  of  the  strips  and  via.  The  current  J 
satisfies 

7s  •  J  -  0  (8) 

and 

n  •  B(J)  -  0  on  S  (9) 

In  (8),  V  •  J  is  the  surface  divergence  of  J.  In  (9),  n  is  a  unit 

8  ~  " 

vector  normal  to  S^,  and  ]}(J)  is  the  magnetic  field  due  to  the  ex¬ 
istence  of  J  in  the  presence  of  the  ground  plane.  Equation  (8)  is  written 
with  the  understanding  that,  at  any  edge  of  that  is  not  at  (z*"±°°),  J 
must  be  tangent  to  that  edge.  If  is  multiply  connected  because  the  via 
has  holes  in  it,  then  the  surface  integral  of  (9)  must  be  zero  over  each 
of  the  holes  [3,  Eq.  (43)].  The  total  current  associated  with  £  is  I. 

In  Fig.  1,  I  flows  from  (z  ■  -<*>)  on  strip  a  through  the  via  to  (z  *■“>)  on 
strip  c.  For  this  reason,  J  must  have  a  z  component  at  the  edges  of 
at  (z“±°°).  Therefore,  can  not  be  tangent  to  the  edges  of  at 
(z«±oo).  Far  from  the  via  on  strip  a,  £  reduces  to  the  current  density 
which  would  exist  on  strip  a  if  both  strip  c  and  the  via  were  absent 
and  if  strip  a  extended  from  (z*-<»)  to  (z  ■  °°) .  Far  from  the  via  on 
strip  c,  J  reduces  to  the  current  density  which  would  exist  on  strip  c 
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if  both  strip  a  and  the  via  were  absent  and  if  strip  c  extended  from 
(2- -oo)  to  (z  ■  ").  If  the  value  of  I  is  specified,  the  conditions 
stated  in  this  paragraph  should  suffice  to  determine  J.. 

In  the  following  manner,  the  integral  over  C  in  (7)  is  transformed 
into  an  Integral  over  S  where  S  is  the  surface  of  the  portion  of  the 
strips  and  via  for  which  (-£  £  z  <_  £c).  Substitution  of  (3)  for  B(J) 
in  (9)  gives 


n  •  V  x  A(J)  -0  on  So,  (10) 

Equation  (10)  is  also  true  on  S  because  S  is  part  of  S^.  Since  (10) 
holds  on  S,  there  is,  according  to  (A-2),  a  scalar  function  f  on  S 
such  that 

A  (J)  -  Vf  on  S  (11) 

—tan  —  s 

where  tan  denotes  the  component  tangent  to  S  and  is  the  surface 
gradient  of  4*  on  S.  Because  of  (11),  expression  (7)  for  becomes 


(V) 


d£ 


Equation  (12)  reduces  to 


to  -  V  -  V 
y  c  a 


(12) 


(13) 


where  V  is  the  value  of  f  at  the  beginning  of  C  at  (z  -  -  t  )  on 
a  a 

strip  a,  and  4*  is  the  value  of  V  at  the  end  of  C  at  (z  ■  on 
strip  c. 

If  l  is  sufficiently  large,  then  A(J)  is  z  directed  on  the 
a 

(z«-£  )  line  on  strip  a,  and  V  V  is  also  z  directed  there  so  that  V 
a  s 


Is  constant  on  the  (z«-£  )  line  on  strip  a.  If  £  is  sufficiently 

a  c 

large,  then  ACJ)  is  z  directed  on  the  (z  ■  ££)  line  on  strip  c,  and 
V  f  is  also  z  directed  there  so  that  4*  is  constant  on  the  (z  *  £  ) 

8  C 

line  on  strip  c. 

Let  J  be  an  electric  current  density  on  S  that  satisfies 


V  •  J  -  0 
s  — 


(14) 


Equation  (14)  is  written  with  the  understanding  that,  at  any  edge  of 

A 

S  that  is  neither  at  (z  =  -£  )  nor  at  (z  ■  £  ),  J  is  tangent  to  that 

a  c  — 

edge.  It  is  assumed  that  J’s  total  current  flowing  from  left  to  right 

A 

in  Fig.  1  is  I.  For  this  reason,  J  must  have  a  z  component  at  the 

edges  of  S  at  (z**-£  )  and  at  (z  *  £  ).  Therefore,  J  cannot  be  tangent 

a  c  — 

to  the  edges  of  S  at  (z*-£  )  and  at  (z  »  £  ).  If  £  and  £  are  suf- 

a  c  a  c 

ficiently  large  so  that  4*  is  constant  on  the  (**”^a)  line  on  strip  a 

and  on  the  (z  •*  £  )  line  on  strip  c,  then  [2,  Eq.  (42)  on  p.  503] 
c 


(4»  -  f  )I  -  f 

c  a  JJ 
S 


7  •  (4M)ds 

s  — 


(15) 


where  4*  is  the  value  of  4*  on  the  (z«-£  )  line  on  strip  a,  and  4* 
a  a  c 

is  the  value  of  4*  on  the  (z  -  ££)  line  on  strip  c.  Substituting 
(4f  -  4»  )  of  (15)  into  (13),  we  obtain 


HI 


7  •  (4'J)ds 

s 


(16) 


which  becomes  [2,  Eq.  (29)  on  p.  502] 


•ill 


(7  4')  •  J  ds  + 


ill 


4*  (7  •  J)ds 


(17) 
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Equations  (11)  and  (14)  reduce  (17)  to 


U»  -  Y  f  A(J)  •  J  ds 


In  contrast  with  (18)  >  the  expression  for  in  terms  of  stored 


energy  is  (1,  Eq.  (4)] 


^  “  Y  jj  AQ.)  *  J  ds 


The  current  J  is  expressed  as 


J  «  J  ,  +  J 

—  — abc  — e 


where 


1J  ,  on  strip  a,  -»  <  z  <  z 
—a  a 

J,  ,  on  the  via,  z  <  z  <  z  (21) 

— b  a  —  —  c 

,  on  strip  c,  <  z  <  00 

In  (21),  J  is  the  current  density  that  J  tends  toward  far  from  the 

3  — 

via  on  strip  a,  and  is  the  current  density  that  J[  tends  toward  far 

from  the  via  on  strip  c.  Since  the  total  current  of  J  is  I,  the  total 

current  of  J  is  I,  and  the  total  current  of  J  is  I.  In  Fig.  1,  I 
—a  — c 

flows  from  left  to  right.  In  (21),  J,  is  any  solenoidal  via  current 

— D 

density  that  is  tangent  to  the  edges  of  S  and  that  provides  for  a  con¬ 
tinuous  flow  of  current  from  J  on  line  a  through  the  via  to  J 

3  t 

on  line  c.  Line  a  is  the  line  at  (z*>z  )  where  the  via  connects  with 

a 

strip  a.  Line  c  is  the  line  at  (z“zc)  where  the  via  connects  with 
strip  c.  Since  the  current  flows  continuously  across  line  a,  no 


11 


electric  charge  can  accumulate  on  line  a  so  that,  on  line  a,  the 

component  of  J,  normal  to  line  a  entering  the  via  la  equal  to  the 
— t> 

z  component  of  J  .  Similarly,  the  component  of  J,  normal  to  line  c 
—a  — b 

leaving  the  via  Is  equal  to  the  z  component  of  J. .  With  J  ,  J.  ,  and  J 

— c  — a  — b  — c 

so  defined,  Is  solenoidal  and  is  tangent  to  any  edge  of  not  at 

(ZB±«>),  The  total  current  of  J  .  is  I.  This  current  flows  continu- 

— abc 

ously  from  (z“-°o)  on  strip  a  through  the  via  to  (z  *  <*>)  on  strip  c. 

Since  both  J  and  J  ,  are  solenoidal,  are  tangent  to  any  edge  of  S  not 
—  — abc  oo 

at  (z*±°°),  and  have  a  total  current  of  I,  the  excess  current  J  in  (20) 

— e 

is  solenoidal,  is  tangent  to  any  edge  of  S^  not  at  (z  ■  ±°°) ,  and  has  no 

total  current  associated  with  it.  Moreover,  J  tends  toward  zero  far 

— e 

from  the  via  on  strip  a  and  far  from  the  via  on  strip  c. 

Substituting  J  .  for  J  and  (20)  for  J  in  (18),  we  obtain 
—abc  —  — 


T  II  * 


A(J  .  )  •  J  ,  ds  + 
—  -abc  —abc 


UN 


-(4}  *  4bcds 


In  the  following  manner,  an  alternative  to  (21)  is  obtained  for  substitution 

Into  the  first  integral  in  (22).  In  (21),  J  is  the  current  density  that 

—a 

would  exist  on  strip  a  if  both  strip  c  and  the  via  were  absent  and  if  strip 

a  extended  from  (z  *  ~°°)  to  (z«°°).  In  (21),  is  the  current  density  that 

would  exist  on  strip  c  if  both  strip  a  and  the  via  were  absent  and  if  strip 

c  extended  from  (z--°°)  to  (z  -  °°) .  Taking  J  and  J  to  extend  from 

a  c 

(z  ■  -°°)  to  (z  ■  °°) ,  we  have 

J  -  J-  +  J+  (23) 

—a  —a  —a 


J  •  J  +  J 
— c  — c  — c 
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where 


J+  -  J 
—a  —a 


J  -  J 
— c  — c 


J+  -  J 
— c  — c 


—  *>  <  Z  <  z 


2  <  z  <  oo 

a  — 


-  CO  <  z  <  Z 

—  C 


z  <  z  < 
c 


(25a) 


(25b) 


(25c) 


(25d) 


Using  (25a)  and  (25d)  and  assuming  that  J.  exists  only  on  the  via*  we 

— b 


can  recast  (21)  as 


J  .  ■  J  +  J.  +  J 

— abc  —a  — b  -c 


Substitution  of  (26)  for  Ja^c  in  the  first  integral  in  (22) 


yields 


(A(Jl+J+)  •  j“  +  A(J"+J.+J+)  •  J.  +  ACT+J,  )  •  J+)ds 
Y  I  JJ  —  — b  — c '  -a  -'—a  -b  -c  -b  -  -a  — b  -c 


+  T  +  T  +  =r  A(J  )  •  J  ,  ds 
a  c  I  —— e  —abc 


where 

T  -iff  A(J~)  •  J~  ds 
a  I  JJ  —  — a  -a 

S 

T  =4  A(J+)  •  J+  ds 

c  I  JJ  '  c /  -c 

S 

Substituting  the  solution  of  (23)  for  J  in  the  argument 
of  A  in  (28),  we  obtain 

T  -  i  ff  A(J  )  •  J~ds  -iff  A(J^)  •  J~  ds 
a  I  JJ  - -a'  -a  I  Jj  -  -a  -a 


In  the  first  Integral  in  (30),  the  domain  of  J  is  extended  to  Include 

A 

(z  ■  z  ) ,  and  J  is  taken  to  be  equal  to  J  at  (z  ■  z  ) .  Now,  the 
first  integral  in  (30)  is  viewed  as  being  over  the  extended  domain  of 
J  rather  than  over  all  of  S.  From  this  viewpoint,  the  first  term  in 
(30)  is  expression  (18)  with  S  replaced  by  the  portion  of  strip  a  for 
which  (-A  <  z  <  z  ) ,  J  replaced  by  J  ,  and  J  replaced  by  J  .  With 
these  replacements,  (18)  is  valid  representation  for  the  flux 
(A  +  z  )  ip  where  \p  ,  the  magnetic  flux  per  unit  length  of  strip  a 
far  from  the  via,  appears  in  (1).  With  its  first  term  so  disposed,  (30) 
reduces  to 

4  *  «a  +  •«>♦.  -  T  H  *<4>  •  4  ds  <31> 

S 

V 

Substituting  the  solution  of  (24)  for  in  the  argument  of 
A  in  (29),  we  obtain 

Tc  "  Y  II  ±<4>  •  4  d"  -  T  II  A<4>  •  4  d8  <32> 

s  s 

In  the  first  integral  in  (32),  the  domain  of  is  extended  to  Include 

(z  ■  z  ),  and  J+  is  taken  to  be  equal  to  J  at  (z  ■  z  ).  Now,  the  first 
C  “C  c 

Integral  in  (32)  is  viewed  as  being  over  the  extended  domain  of  rather 
than  over  all  of  S.  From  this  viewpoint,  the  first  term  in  (32)  is  ex¬ 
pression  (18)  with  S  replaced  by  the  portion  of  strip  c  for  which 

A  ^ 

(z  <  z  <  A  ),  J  replaced  by  J  ,  and  J  replaced  by  J  .  With  these  re- 
c  —  —  c  —  — c  —  C 

placements,  (18)  is  a  valid  representation  for  the  flux  (A£  -  Zj.)^ 

where  ijs  ,  the  magnetic  flux  per  unit  length  of  strip  c  far  from  the  via, 
c 
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appears  In  (1).  With  Its  first  term  so  disposed,  (32)  reduces  to 


Tc-  II  A<£>  -4ds 

s 


Substituting  (31)  and  (33)  Into  (27)  and  then  using  (26),  we 


obtain 


-  T  I  [A(J+-J++J  +J.  )  •  J~  +  A(J“-J"+J  +J.  )  •  J+  +  A(J"+J++J  «_  )  •  JL  )ds 
1  iJ - c-a-e-b  -a - a-c-e-b  -c - a-c-e-b  -b 


+  (A  +  *)’/'_  +  (i  ~  z  )ij> 
a  a  a  c  c  c 


When  (34)  Is  substituted  Into  (1),  the  flux  (i  if/  +  Z  ip  )  cancels  and  we 

A  A  C  C 

are  left  with 

L  -  11m  1 1  (A(J+- J++J  +J.  )  •  J>A(J~-J>J  +J.)«/+*(jV+J  +J.  )  -  JL)da 

«  H  -ko  i2  J  J  -a  -e  -b  -*  -  -a  -c  -e  -V  -c  -'-a  -a  -e  -b' 

a  S 


z  if/  -  i  if) 

a  a  c  c 

I 


where  S  Is  the  surface  of  the  portion  of  the  strips  and  via  for  which 

(■^  £  z  £  A  )  • 
a  —  —  c 

If  £  increases  by  the  amount  A£  ,  then  the  quantity  in  brackets  on 

A  A 

the  right-hand  side  of  (35)  increases  by 


pr  ||  *<4  -  4  +  4  +  4>  -4,i5 


where  is  the  portion  of  strip  a  for  which  (-A#  £  z  £  “  £g)' 

If  8,  Is  large  and  if  Af.  is  positive,  then  the  vector  potential  A 
a  a 
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in  (36)  is  small  on  because  its  current  source  is  distant  from 

**A?a'  Moreover,  it  is  conjectured  that,  if  P-a  is  sufficiently  large, 

then  expression  (36)  is  negligibly  small  for  (0  <  A£  -^  “»). 

a 

If  increases  by  the  amount  AP.^,  then  the  quantity  in  brackets  on 
the  right-hand  side  of  (35)  Increases  by 


i-  ff  A(.r  - 

2  JJ  --a 


A  ?.c 


4+  4> 


.I+ds 


(37) 


where  S,„  is  the  portion  of  strip  c  for  which  (fc  •'  z  P  +  A?  ) . 

A  tc  r  -  —  c  c 

If  9  is  large  and  if  AP  is  positive,  then  the  vector  potential  A 
c  c  — 

in  (37)  is  small  on  ^  because  its  current  source  is  distant  from 
Moreover,  it  is  conjectured  that,  if  9.  is  sufficiently  large,  then 
expression  (37)  is  negligibly  small  for  (0  Alt  '  •") . 

If  the  conjectures  in  the  previous  two  paragraphs  are  true,  then 
the  limit  in  (35)  exists.  Tn  this  case,  a  good  approximation  to  Lg  is  given 
by 


75 II  <*<4-4+ 


.T  +.f)  •  .I'+Ad'-J'+J  +.1.)  •  J++A(J"+J++J  +JL)  •  Ju)ds 
-e  ~b  -a  --  -a  -c  -e  — b  — c - a  -c  — e  — b  — b 


t.  •!>  -  t  if 

a  a  c  c 
*  - 


(38) 


where  S  is  the  surface  of  the  portion  of  the  strips  and  via  for  which 

♦ 

(— V,  z  *  l  )  where  9  and  V.  are  sufficiently  large.  In  (38),  J~  and 
a  —  —  c  ac  -a 

J~  are  given  by  (25).  Furthermore,  .1  ,  J  ,  -1  ,  and  J,  are  described  in 
— c  —a  — c  - e  — 0 

the  paragraph  that  contains  (21).  In  (38),  MJ^-J^+J+J^)  is  the  magnetic 
vector  potential  due  to  the  existence  of  +J  +J.  )  in  the  presence 

C  8  "“G  O 


of  the  ground  plane.  The  remaining  A's  in  (38)  are  similarly  defined. 


Moreover,  ip^  is  the  magnetic  flux  per  unit  length  of  strip  a  far  from 

the  via,  and  ip is  the  magnetic  flux  per  unit  length  of  strip  c  far  from 

the  via.  Otherwise  stated,  if^/I  is  the  inductance  per  unit  length  of 

strip  a  far  from  the  via,  and  ijWl  is  the  inductance  per  unit  length  of 

strip  c  far  from  the  via.  In  Section  III,  the  fluxes  <p  and  ip  will 

a  c 

emerge  as  by-products  of  the  solutions  for  and  .J^.  The  right-hand 

side  of  (38)  does  not  depend  on  I  because  the  A's,  the  J.'s,  and  the  i|»'s 
in  (38)  are  proportional  to  1. 


III.  INTEGRAL  EQUATION  FORMULATION 

The  objective  of  Section  HI  is  to  derive  integral  equations  for  J  , 

— a 

.1  ,  and  J  and  to  obtain  expressions  for  ip  and  ip  . 

— c  —e  r  a  c 

By  definition,  J  is  the  current  density  that  J  tends  toward  far  from 

—a  — 

the  via  on  strip  a.  Equivalently,  is  what  J  would  be  if  both  strip  c 

and  the  via  were  absent  and  if  strip  a  extended  from  (z  -  -  ntj  to  (z  ■  °°) , 

Hence,  the  problem  of  determining  J  is  a  two-dimensional  one.  The  situ- 

d 

ation  is  shown  in  Fig.  j.  Since  J's  total  current  flowing  from  left  to 
right  in  Fig.  1  is  I,  J  's  total  z-directed  current  is  l. 

— 3 

Because  d  satisfies  (8)  and  (9),  must  satisfy 


and 


v  •  J  *0  (39) 

s  —a 


u 


y 


B(J  ) »  0  on  strip  a 
—a 


(40) 


In  (40),  u  is  the  unit  vector  in  the  y  direction,  and  B(J  )  is  the 

***Y  "fl 


0 


x 


x 

o  a 


An  infinitely  long  perfectly  conducting  strip  of  width 
w  at  height  h  above  a  perfectly  conducting  ground  plane. 

fl  3 

On  the  strip,  there  is  a  magnetostatic  current  density 
whose  total  z-directed  current  is  1. 
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magnetic  field  due  to  the  existence  of  J  in  the  presence*of  the  ground 

8 

plane.  Equation  (39)  is  written  with  the  understanding  that,  at  any  one 

of  the  two  lateral  edges  of  strip  a,  ±  must  be  tangent  to  that  edge. 

The  magnetic  field  B(J  )  that  appears  in  (40)  is  expressed  as 
1  a 


B(J  )  =  V  x  A(J  )  (41) 

where  A(J  )  is  the  magnetic  vector  potential  due  to  the  existence  of  J 
- a  8 

in  the  presence  of  the  ground  plane.  If  J  was  an  arbitrarily  directed 

3 

current  above  the  ground  plane,  then 


A(J  ) 
- a 


_LL 

4ir 


II 


strip  a 


4<I')  .  v,y(i'>  - 
ir^rr  ir^n 


u  J  (r*) 
— z  az  — 


)  ds ' 


(42) 


where  (u  ,  u  ,  u  )  are  the  unit  vectors  in  the  (x,y,z)  directions,  and 
— x  — y  —  z 

(J  ,  J  ,  J  )  are  the  (x,y,z)  components  of  J  .  In  (42),  ds'  is  the  dif¬ 
ax  ay  az  —a 

ferential  element  of  area  at  r'  on  strip  a,  r"  is  the  image  of  jr *  about  the 
ground  plane,  and  jr  is  the  point  at  which  ACJ^)  is  evaluated. 

Because  the  problem  of  determining  J  is  two-dimensional,  J  is  inde- 
pendent  of  z  so  that 


J 

—a 


J 

ax 


<x)ux 


+  J  (x)u 
az  -z 


(43) 


Now,  (39)  expands  to  [2,  Eq.  (18)  on  p.  501] 


DJ  (x)  3J_  (x) 

AX  ^  A  2 


dX 


(44) 


Since  J  does  not  depend  on  z,  (44)  implies  that  J  does  not  depend  on 
az  ax 

x.  Moreover,  J  is  zero  at  (x  -  x  )  and  at  (x  -  x  +w  )  because  J  was 

aX  3  a  a  ■ 


constrained  to  be  tangent  to  the  lateral  edges  of  strip  a.  Therefore* 


is  zero  everywhere  on  strip  a.  Hence,  (43)  simplifies  to 


J  ■  J  (x)u 
—a  az  —t 


(45) 


Equation  (45)  reduces  (42)  to 


A(J  )  -  A  (J  )u 
- a  z  —a  — z 


(46) 


where 


x  +w 
a  a 


V4 >  -  &  J  j  d2'/ 


/(x-x * ) 2+(y-ha) 2+(z-z ’ ) 2  / (x-x,)2+(y+h^^-(z-z’/ 


(47) 


In  (47),  (x.y.z)  are  the  rectangular  coordinates  of  the  point  at  which 
A  (J  )  is  evaluated. 


If  the  Integral  with  respect  to  z'  in  (47)  is  called  I  then  [4, 


Formula  200.01] 


1  -  lim 

a-K» 


g-x» 


+  V Y2  +  (x-x’)2+  (y-h  )2 


y-e 


In  ( 


Y  +  /y2+(x-x')2  +  (y  +  h  )2  J 


(48) 


Y”  -a 


The  logarithm  in  (48)  vanishes  at  its  upper  limit  so  that  (48)  reduces 

to  _ 

- a+  /a2  +  (x-x*)2  +  (y+h  )2 

(49) 


I  ..  “lim  In  ( 
zl  a- 


=) 


-  a  +  ^a2  +  (x-x*)2  +  (y-h  )2 


Expression  (49)  is  rewritten  as 
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(x-x')2  +  (y+h  )2  a  +  Ju2  +  (x-x')2  +  (y-h  )2 

■  llm  ln  ( - 5 - ^-sr)  ( - - - 9—: 

cr**5  (x-x')  +  (y-h  )  f~2  2  5 

L  y  a;  a  +  yoi  +  (x-x')2  +  (y+h  )2 

fl 


)  (50) 


The  function  of  a  in  parentheses  in  (50)  approaches  unity  so  that  (50) 


becomes 


(x-x')2  + (y+h  )2 
lzi  =  ln  ( - t - ^ 

(x-x' )  +  (y-h  y 
a 

Substituting  expression  (51)  for  the  integral  with  respect  to  z'  in 


(47),  we  obtain 


x  +w 
a  ri 


A  (J  ) 
z  —a 


_U_ 

4-n 


3  (  1  (x-x')2  +  (y-fh  )2 

Jaz(x,)  ln  ( - 2 - V)  dx'  (52) 

J  (x-x')2+  (y-h  )2 

Xa  3 


In  view  of  (46),  substitution  of  (41)  into  (40)  gives 


u  •  V  *  (A  (J)u)  *  0  on  strip  a 
y  z  a  z 

Equation  (53)  reduces  to  (2,  Eq .  (51)  on  p.  490] 


M  (J  ) 
z  —a 


0  on  strip  a 


It  is  evident  from  (52)  and  (54)  that  Az(Ja)  is  constant  on  strip  a. 
Hence,  using  (52),  we  have 


x  +w 
a  a 


(x-x’)2  +  (2h  )2 


J  (x  )  ln  ( - r - )  dx'  =  u,  x  <  x  x  +w  (55) 

az  .,.2  a  —  —  a  a 

(x-x') 


where  a  is  a  constant  that  is  determined  by  requiring  J  's  total  z- 

az 


directed  current  to  be  I. 


The  flux  \p  is  given,  in  analogy  with  (7),  by 
a 


)  •  dJ 


where  the  line  integral  in  (62)  is  over  a  unit  length  of  strip  a  in 
the  z  direction.  Substituting  (46)  into  (62),  we  obtain 


»  a  (j  ) 

a  z  — a 


where  A  (J  )  is  evaluated  on  strip  a.  Expression  (52)  gives 
z  — a 


X  +w 
a  a 


A  (J  )  -  J  (x')  In  ( 

z  —a  4^i  az 


(x-x')2  +  (2h  )2 

_ S—i 

2  ‘ 

(x-x'T 


-)dx'  on  strip  a 


Substituting  (61)  into  (64)  and  using  (57)  and  (63),  we  obtain 


x  +w 
a  a 


P 

1  (x')dx' 

az 


where  J  (x')  satisfies  (57). 
az 

The  inductance  per  unit  length  of  strip  a  is  called  L  and 

8 

is  defined  by 


La-f 


Substituting  (65)  into  (66),  we  obtain 


^a  x  +v 

fa  a 


J  (x ' )  dx  ' 
az 


where  J  (x')  satisfies  (57) 
az 
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The  capacitance  per  unit  length  of  strip  a  in  [5]  is  called 

C  and  is  defined  by 
a  1 


"a  V 


where 


x  +w 
a  a 


q  (x  ' ) dx ' 
a 


In  (69),  q  (x')  satisfies  (5,  Eq.  (9)].  Substitution  of  (69)  into 

fl 


(68)  gives 


x  +w 
a  a 

J  q, (*•>*■• 

X 

a 


where  q  (x ' )  satisfies  [5,  Eq.  (9)].  If 
3 


where  e  is  the  permittivity  above  the  ground  plane  in  [5],  then  com¬ 
parison  of  [5,  Eq.  (9)]  with  (57)  gives 


9  (x')  -  J  (x')  (72 

»  mZ 

Because  of  (71)  and  (72),  the  product  of  the  right-hand  sides  of  (67) 
and  (70)  is  pe  so  that 


L  C  -  pe 
a  a 


Relationship  (73)  is  well-known  for  a  two  conductor  lossless  trans¬ 
mission  line  [6,  Eq.  (10b)  on  p.  123]. 

By  definition,  J  is  the  current  density  that  J  tends  toward 
— c 

far  from  the  via  on  strip  c.  Equivalently,  J  is  what  J  would  be  if 
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both  strip  a  and  the  via  were  absent  and  if  strip  c  extended  from 
(z»-oo)  to  (z  =*  °°) .  Hence,  the  problem  of  determining  is  a  two- 
dimensional  one.  The  result  of  a  development  similar  to  (39)-(61) 
is  that 


where 


J  =  J  (x)u 
— c  cz  — z 


I  J  (x) 

j  (x) - - - “ - 

CZ  X  +w 

c  r  c 


J  (x ' )dx ' 
cz 


where  J  (x1)  satisfies 
cz 


x  +w 
c  c 


J  (x')  In  (- 
cz 


(x-x’)2  +  (2h  )2 


- 5 - )  dx'  «  1,  x  _<  x  <  x  +w  (76) 

(x-x!) 


The  analogue  of  (65)  is 


x  +w 
c  c 


J  (x')dx' 
cz 


where  J  (xf)  satisfies  (76). 
cz 

Now,  the  objective  is  to  derive  an  integral  equation  for  J^. 

In  (20),  is  solenoidal,  is  tangent  to  all  the  edges  of  except 
those  at  (z=±OT),  and  its  total  current  flowing  from  left  to  right 
in  Fig.  1  is  I.  Moreover,  is  solenoidal,  is  tangent  to  all  the 

edges  of  S  except  those  at  (z=±°°),  and  its  total  current  is  I. 
Therefore,  J  is  solenoidal,  is  tangent  to  all  the  edges  of  with  the 
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possible  exception  of  those  at  (z=±°°),  and  has  no  total  current  associ¬ 
ated  with  it.  Since  tends  toward  zero  far  from  the  via  on  strip  a 

and  far  from  the  via  on  strip  c,  J  is  zero  on  the  edges  of  S  that  occur 

— C  oo 

at  (z“±°°).  Hence,  J  is  tangent  to  all  the  edges  of  S  .  The  solenoidal 

— Tg  00 

property  of  is  expressed  as 

Vs  *  4  -  0  <*«) 

Substitution  of  (20)  into  (9)  gives 

n  •  B(J  )  -  -  n  •  B(J  .  )  on  S  (79) 

— - e  — - abc  00 

where  JiOJ^)  is  the  magnetic  field  due  to  the  existence  of  in  the 
presence  of  the  ground  plane,  and  B(Jabc)  is  the  magnetic  field  due  to 
the  existence  of  in  the  presence  of  the  ground  plane. 

Because  of  (78),  there  exists  a  scalar  function  u(r)  such  that  [3, 
Eq.  (B-l)l 

J  (r)  *  n  x  V  u(r)  (80) 

— e  —  —  s  — 

where  V  is  the  surface  gradient  on  S^.  Since  J  is  tangent  to  each 
edge  of  S^,  (80)  implies  that  u  is  constant  on  each  edge  of  S^.  If  the 

via  is  a  simple  strip,  then  S^  has,  in  the  domain  of  finite  z,  only  two 

lateral  edges.  In  this  case  (80)  predicts  that  the  value  of  u  will  be 

the  same  on  both  of  these  edges  because  J  has  no  total  current  associ- 

— e 

ated  with  it.  In  general,  however,  the  value  of  u  on  one  edge  of  is 
not  necessarily  the  same  as  that  on  any  other  edge  of  S^.  Only  deriva¬ 
tives  of  u  appear  in  (80).  Hence,  if  (80)  is  substituted  into  (79), 
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then  (79)  can  only  determine  u  to  within  an  arbitrary  additive  constant. 
To  make  u  unique,  we  specify  (u  =  0)  at  an  arbitrary  point  on  S  .  If 
the  via  is  a  simple  strip,  it  is  natural  to  specify  (u  *  0)  at  a  point 
on  one  of  the  two  lateral  edges  of  because  this  specification  will 
render  (u  «  0)  on  the  entirety  of  both  lateral  edges  of  S^.  In  summary, 
the  boundary  conditions  on  u  in  (80)  are  that  u  is  constant  on  each 
edge  of  and  that  (u  =  0)  at  some  point  on  S^. 

In  (79),  B(J  )  is  expressed  as 


B (J  )  =  V  x  A(J  ) 


where 

re  J  (r')  u  J  (r*)  -  u  J  (r’)  -  u  j  (r') 

.  ,T  W  ,-e  —  ,  -y  ey  -  -x  ex  -  -zJez  -  ,  ,  , 

"47  tT^T  +  - rr-"r7rT  - 1  d8  (82) 


where  (J  ,  J  ,  J  )  are  the  (x,y,z)  components  of  J  .  In  (82),  jr' 

is  the  point  at  which  the  differential  element  of  area  da'  is  located, 

r"  is  the  image  of  jr'  about  the  ground  plane,  and  r_  is  the  point  at 

which  A(J  )  is  evaluated.  Equations  (81)  and  (82)  are  also  valid  with 
- e 

J  replaced  by  J  ,  .  With  li  given  by  (81),  (79)  becomes 

G  9.  DC 


n  •  V  xA(J  )  =  -  n  •  V  x  A(J  ,  )  on  S 
—  —  — e  —  —  —a  be  a 


where  A  is  given  by  (82),  and  is  given  by  (80)  in  which  u  is  an 
unknown  scalar  that  is  constant  on  each  edge  of  and  whose  value  is 
zero  at  some  point  on  SM.  If  is  multiply  connected  because  the  via 
has  holes  in  it,  then  (83)  must  be  accompanied  by  the  auxiliary  condition 


that  the  surface  integral  of  the  left-hand  side  of  (83)  over  any  one  of 
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the  holes  must  equal  the  surface  integral  oi  the  right-hand  side  of 
(83)  over  that  hole. 

Equation  (83)  represents  a  surface  integral  equation  for  u.  If 
(83)  can  be  solved  for  u,  then  will  he  given  hy  (80).  Rigorously, 
the  domains  of  both  u  and  the  point  at  which  both  sides  of  (83)  are 
observed  extend  to  (z*-“)  in  negative  z  direction  a  in)  to  (z  =  •»’)  in 
the  positive  z  direction.  However,  as  |z|  glows,  V  u  approaches 
zero,  and  both  sides  of  (83)  approach  zero.  Since  'Mi  approaches 
zero  as  |z|  grows,  the  point  at  which  (u  -  0)  can  be  chosen  so  that 
u  approaches  zero  as  |z|  grows.  In  Section  IV,  a  numerical  solution 
for  u  is  obtained  by  truncat ing  both  the  unknown  u  and  its  integral 
equation  at  a  finite  negative  value  of  z  and  at  a  finite  positive  value 
of  z . 


IV.  SOLUTION  BY  THE  METHOD  Of  MOMENTS 

In  Section  IV,  the  method  of  moments  is  used  to  numerically  solve 

the  integral  equations  (57),  (76),  and  (83)  for  J  ,  and  u»  re- 

spectively.  After  these  equations  have  been  solved,  can  be  obtained 

by  substituting  (61)  into  (45),  can  be  obtained  by  substituting  (75) 

into  (74),  and  J  can  be  obtained  by  simply  substituting  u  into  (80). 

-e 

The  quantities  I  ,  J+,  J  ,  J  +  ,  <j'  ,  and  ip  in  expression  (38)  for  the 

—a  —a  — e  -c  a  c 

excess  inductance  I.  will  then  he  given  by  (25a),  (25b),  (25c),  (25d), 
e 

(65),  and  (77),  respectively. 

The  moment  solution  for  .1  is  constructed  by  expanding  J 

az  az 

as  a  linear  combination  of  pulse  functions  and  by  point  matching  (57) 


m 
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at  the  centers  of  the  doaaina  of  the  pulses.  The  sonant  solution  for 

A 

J  is  constructed  in  the  seat  manner.  The  moment  solution  for  u  is 
cs 

constructed  by  truncating  at  (s  ■  -  on  the  negative  s  axis  and 
at  (z  •  &c)  on  the  positive  z  axis,  calling  the  truncated  surface  S, 
modeling  S  by  triangular  patches  [ 7] •  calling  the  resulting  triangu¬ 
lated  surface  S^,  expanding  u  as  a  linear  combination  of  pyramid 
functions  on  ST»  and  testing  (83)  with  the  pyramid  functions.  Such 
testing  consists  of  multiplying  (83)  by  each  of  the  pyramid  functions 
successively  and  integrating  over  S^.  In  the  resulting  simultaneous 
equations,  the  derivatives  on  A  are  transferred  to  the  pyramid  testing 
functions,  and  the  Integrals  over  the  pyramid  testing  functions  are 
approximated  by  sampling  A  at  the  centroids  of  the  triangular  patches. 

To  construct  the  moment  solution  of  (57)  for  J^.  we  approximate 


J 

az 


by 


I 

where 


(84) 


X.J  -  X.  +  J-1.2 . V1 


(85) 


A  ■  w  /N  (86) 

a  a  a 

In  (84).  P  (x)  is  the  pulse  function  defined  by 
a 

^l,  0  <  x  <  Ag 

P  (x)  "  i  (87) 

8  I 

■  0,  otherwise 

Not  used  in  (84),  the  quantity  x  „  ..  of  (85)  is  the  x  coordinate 
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of  the  farther  edge  of  strip  a.  This  quantity  is  used  later. 

Next,  N  equations  involving  the  unknown  constants  {i  .}  are 
a  aj 

obtained  by  substituting  (84)  into  (57)  and  enforcing  (57)  at  the 
centers  of  the  domains  of  the  pulses  in  (84).  These  equations  are 
written  in  matrix  notation  as 


P  I  =  V 
a  a  a 


(88) 


where  I  and  V  are  column  vectors  The  jth  element  of  I  is  I  . ,  and 
a  a  a  aj 

each  element  of  V  is  unity.  In  (88),  P  is  a  square  matrix.  Its  ijth 
3  3 

element  is  given  by 


aij 


a,j+l 

j 


aj 


(x+.  -  x’)2  +  (2h  )2  '  1=1,2 . N- 

ln  ( — — - - - g-A-jdx'  . 

(xai  "  X<)  /  j=l  ,2 . N 


(89) 


where 


Xai  =  Xa  +  (i  ‘  I)Aa 


(90) 


Expression  (89)  for  P  becomes  [5,  Eq .  ( B— 8 > ] 

aij 


aij 


2h  2  , 

x  In  (1  +  (— )  )  +  4h  tan  (^-) 

x  3  In 

a 


(lJ-il+2  )\ 


(91) 


The  moment  solution  of  (57)  for  J  is  completed  by  solving  (88)  for 

3Z 

I  and  substituting  the  elements  of  I  into  expression  (84)  for  J  . 

a  3  3Z 

Similar  to  the  moment  solution  of  (57)  for  J  ,  the  moment  solution 


of  (76)  for  .1  is  given  by 
cz 


N 


^  ,  ,x  V  I  .P  (*'  -  X  .) 

.1  (x  )  =  l  C]  c  cj 

cz  j=l 


(92) 
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where 


x  .  »  x  +  (j-l)A  ,  J*l,2 . N  +1 

C  J  c  c  c 


L  *  w  /N 
c  c  c 


PcOO  = 


1,  0  <  x  A 


otherwise 


Not  used  in  (92),  the  quantity  x  .  of  (93)  is  the  x  coordinate  of 

c,N  +1 

c 

the  farther  edge  of  strip  c.  In  (92),  T  .  is  the  jth  element  of  the 

J 

-f 

column  vector  I  that  satisfies 

c 


PI  *=  V 
c  c  c 


Here,  V  is  a  column  vector  of  N  elements.  Each  element  of  is 
c  c  c 

unity.  In  (96),  Pc  is  a  square  matrix  of  order  .  In  analogy  with 
(91),  the  ijtli  element  of  P^  is  given  by 


2h  , 

P  ,  -  x  In  (1  +  (— £)  )  +  4h  tan_i(^-) 

cij  x  c  zhc 


(  1 1~i !  +!>*. 


(H-il  -f)Ac 


Before  constructing  the  moment  solution  of  (83)  for  u,  we  will 

show  that  the  solution  u  to  (83)  is  proportional  to  I.  Certainly, 

u  will  be  proportional  to  I  if  on  the  right-hand  side  of  (83)  is 

proportional  to  I.  According  to  (21),  Ja^c  consists  of  the  part  of  J_a 

for  which  (-°°  <  z  <z  ),  J.  on  the  via,  and  the  part  of  J  for  which 

a  — b  — c 

(z  <  z  <  °°) .  Because  of  (45)  and  (61),  J  is  proportional  to  I.  Be- 
c  —a 

cause  of  (74)  and  (75),  J  is  proportional  to  I.  Since  J,  is  any 

— C  — D 

solenoidal  via  current  density  that  is  tangent  to  the  edges  of  S  and 
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that  provides  for  a  continuous  flow  of  current  from  Jg  on  strip  a 

to  on  strip  c,  can  be  chosen  to  be  proportional  to  I.  With 

J.  so  chosen,  each  of  the  currents  J  ,  J,  ,  and  J  is  proportional 
— b  —a  — b  — c 

to  I.  Therefore,  J^c  is  proportional  to  I.  Hence,  u  is  propor¬ 
tional  to  I. 

Construction  of  the  moment  solution  of  (83)  for  u  is  begun 


by  expanding  u  as 


u(r)  «  I  [  I.u  (r) 
j  =  l  J  J 


where  (u^ }  are  expansion  functions  and  {i^}  are  unknown  constants  to 

be  determined.  Because  u  was  shown  to  be  proportional  to  I  in  the 

previous  paragraph,  none  of  the  unknowns  { I j }  will  depend  on  I.  Before 

defining  the  expansion  functions  {u.},  we  truncate  strip  a  at  ) 

J  a 

far  from  the  via  and  truncate  strip  c  at  (z  »  9..  )  far  from  the  via.  The 

c 

surface  of  the  via  and  the  truncated  strips  is  called  S  and  is  modeled 
by  triangular  patches  [7],  The  resulting  triangular  patch  surface  is 
called  S,j,.  The  vertices  of  the  triangles  are  called  the  nodes  of  S^,. 
The  triangles  are  chosen  so  that  none  of  them  straddles  strip  a  and 
the  via,  and  none  of  them  straddles  strip  c  and  the  via.  Each  triangle 
is  either  entirely  on  strip  a,  or  entirely  on  the  via,  or  entirely  on 
strip  c. 

If  the  via  is  a  simple  strip,  then  there  is  a  one-to-one  cor¬ 
respondence  between  the  expansion  functions  and  the  interior  nodes  of 


ST»  An  interior  node  of  ST  is  a  node  of  not  on  any  edge  of  S^. 
The  expansion  function  u^  is  associated  with  the  jth  interior  node 


Consider  each  triangle  of  which  one  vertex  is  the  jth  interior 


of  S_. 

T 

node  of  S^.  Called  (T^,  n=l,2, . . . ,Nj ) ,  these  triangles  encircle  the 
jth  interior  node  of  ST*  The  set  of  these  triangles  is  the  domain 
of  Uj .  On  Tjn>  Uj  decreases  linearly  from  unity  at  the  jth  interior 
node  of  to  zero  at  the  side  of  T^n  opposite  the  jth  interior  node 
of  S^,.  More  precisely, 


Oj(£)  -  Cjn(r) 


on  T.  ,  n=l,2 . N, 

jo  j 


(99) 


where  f  .  is  the  area  coordinate  associated  with  the  vertex  of  T, 
jn  jn 

that  is  the  jth  interior  node  of  S_.  At  the  point  r  on  T.  ,  £.  is 

T  jn  jn 

defined  by 


AT. 

C ,  (r )  = 

Jn  -  Ajn 


(100) 


where  A,  is  the  area  of  T,  ,  and  A.  is  the  area  T.  would  have  if  its 
jn  jn  jn  jn 

vertex  which  is  the  jth  interior  node  of  were  replaced  by  the  point  r. 

As  defined  by  (99),  the  expansion  function  is  unity  at  the  jth  interior 

node  of  S^,  is  zero  on  the  edges  of  S^,,  and  is  continuous  everywhere. 

Furthermore,  V  u.  is  continuous  everywhere  except  on  the  sides  of  the 
s  j 

triangles  that  encircle  the  jth  interior  node  of  S  . 

If  J,  is  the  expansion  function  for  J  associated  with  u.,  then 
-j  — e  j 

according  to  (80), 


ij  (r)  =  n  x  Vgu j  (jr) 


(101) 


Substituting  (99)  into  (101),  we  obtain 


l. 

-In 


-i  2A 


jn 


on  T,  ,  n»l ,2, . . . ,N 
jn  J 


(102) 


where  H.  is  a  vector  whose  length  is  the  length  of  the  side  of  T 
-jn 


jn 


opposite  the  jth  interior  node  of  and  whose  direction  is  parallel 

to  this  side.  The  sense  of  £.  is  such  that  J.  encircles  n  in  the 

~jn  -j 

left-hand  sense.  A  left-handed  screw  whose  axis  is  parallel  to  ti  would 

advance  in  the  direction  of  n  when  rotated  in  the  direction  of  JL  . 

—  —  jn 

The  current  Jj  flows  continuously  around  the  jth  interior  node  of  S^,  in 

the  sense  that,  on  each  triangle  side  that  fans  out  from  this  node,  the 

component  of  current  normal  to  this  side  is  continuous.  Here,  the  word 

continuous  is  used  loosely  because,  if  S^,  bends  sharply  at  the  side  of 

a  triangle,  the  component  of  current  normal  to  this  side  must  change 

direction  suddenly  in  order  to  remain  on  ST  .  Obviously,  J  is  tangent 

L  j 

to  all  triangle  sides  that  are  edges  of  ST»  As  defined  by  (102),  J. 

is  J  /2  where  J  is  shown  in  [8,  Fig.  2.2]. 

L  L 

If  the  via  is  not  a  simple  strip,  S^,  may  contain  junctions  of 
surfaces.  Each  junction  of  surfaces  is  a  chain  of  straight  line  seg¬ 
ments  drawn  between  nodes.  This  chain  of  straight  line  segments  is 
called  a  junction  line,  and  these  nodes  are  called  junction  nodes. 


If  some  electrical  contacts  are  ignored  along  junction  lines. 


then  the  junctions  will  disappear,  and  S^,  will  separate  into  several 

isolated  surfaces.  Let  be  a  typical  one  of  these  surfaces.  If 

S.  is  closed,  all  the  nodes  of  S,  are  interior  nodes,  and  an  expansion 
K  K 

function  is  associated  with  each  of  these  nodes  except  one.  If  is 
open,  an  expansion  function  is  associated  with  each  interior  node  of  S^. 
Just  after  (83),  it  was  stated  that  the  value  of  u  is  zero  at  some 
point  on  the  surface.  With  the  expansion  functions  chosen  previously 
in  this  paragraph,  u  is  zero  at  the  edges  of  if  is  open.  If 

is  closed,  u  is  zero  at  the  interior  node  which  does  not  have  an 


34 


expansion  function  associated  with  it.  If  is  open  and  multiply 
connected  because  it  has  holes  in  it,  then  has  more  than  one  edge. 
In  this  case,  additional  expansion  functions  are  needed  because, 
although  it  is  proper  to  have  (u=0)  on  one  edge,  it  is  too 
restrictive  to  have  (u=0)  on  all  edges.  Currents  can  circulate  about 
the  holes.  None  of  these  circulating  currents  can  be  expressed  as  a 
linear  combination  of  the  expansion  functions  associated  with  the  in¬ 
terior  nodes  of  S,  . 

k 

An  additional  expansion  function  must  be  associated  with  each 
hole  of  S^.  A  typical  hole  of  has  a  closed  contour  which  is  a  chain 
of  straight  line  segments  drawn  between  nodes.  This  chain  of  straight 
line  segments  is  called  the  contour  C.^.  The  expansion  function 
associated  with  the  hole  whose  contour  is  Sole  is  called  u^0^e«  The 
domain  of  uho^e  consists  of  all  the  triangles  that  are  attached  to 
Sole '  A  trian8le  is  attached  to  Chole  if  at  least  one  of  its  ver¬ 
tices  is  a  node  of  C^ole.  By  definition,  u^ole  unity  on  clloje* 

The  nodes  of  C.  ,  are  called  C,  .  nodes.  On  each  triangle  of  which 
hole  hole 

exactly  one  vertex  is  a  node,  u^Q^e  decreases  linearly  from  unity 

at  the  C,  .  node  to  zero  at  the  side  opposite  this  node.  On  each 
triangle  of  which  exactly  two  vertices  are  nodes,  u^0^e  decreases 

linearly  from  unity  on  the  line  segment  which  connects  these  nodes  to 
zero  at  the  remaining  vertex.  Finally,  u^o^e  is  unity  everywhere  on 
each  triangle  of  which  all  three  vertices  are  nodes.  The  expan¬ 

sion  function  for  associated  with  u^o^e  is  called  J^o^e  an<*  is  given, 
according  to  (80),  by 
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J,  (r)  =  n  x  V  u  (r) 
—hole  —  —  s  hole  — 


(103) 


If  the  electrical  contacts  that  were  previously  ignored  along 

junction  lines  are  restored,  loops  about  which  current  can  circulate 

are  formed.  Around  each  loop,  a  closed  contour  is  chosen.  For  each  of 

these  contours,  an  additional  expansion  function  is  needed.  Let  be 

a  typical  one  of  these  contours  and  let  be  the  expansion  function 

associated  with  it.  Similar  to  the  expansion  function  defined  in 

the  previous  paragraph,  u.,  is  unity  on  C  ,  and  u,  decreases 

linearly  to  zero  at  the  string  of  adjacent  nodes  that  run  along  one  side 

of  The  expansion  function  u^o^e  was  defined  only  to  one  side  of 

C,  ,  because  C,  ,  was  a  boundary  of  S,  .  However,  C,  is  not  always 
hole  hole  k  loop 

a  boundary  of  S„  so  that  u,  must  be  defined  to  both  sides  of  C,  .  To 
T  loop  loop 

one  side  of  cioop *  ui00p  was  required  to  decrease  linearly  to  zero  at  the 

adjacent  nodes.  To  the  other  side  of  ,  we  want  uioop  to  suddenly 

drop  to  zero.  Unfortunately,  the  ensuing  discontinuity  in  u^oop  will  give 

rise  to  an  impulse  in  7  u,  .  The  expansion  function  for  J  associated 
r  s  loop  — e 

with  u,  is  called  J,  .  According  to  (80),  J..  is  given  by 
loop  —loop  —loop 


J.  (r)  =  n  x  u1  (r) 
—loop  —  —  s  loop  — 


(104) 


However,  not  wanting  any  impulse  in  we  alter  £^00p  80  that 

—loop  is  given’  not  strictly  by  (104),  but  by  (104)  with  Vg  uioop^ 
stripped  of  its  impulse. 

When  the  previously  ignored  electrical  contacts  are  restored  along 
junction  lines,  it  may  be  necessary  to  associate  additional  expansion 
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functions  with  some  of  the  junction  nodes.  A  typical  expansion  function 

associated  with  a  junction  node  is  called  u.  .  Now,  u.  is  either 

junc  June 

an  expansion  function  that  was  associated  with  the  junction  node  before 

the  electrical  contacts  were  restored  or  an  expansion  function  that  must 

be  appended  when  the  electrical  contacts  are  restored.  The  domain  of 

Ujunc  completely  surrounds  the  junction  node,  part  of  the  domain  being 

on  one  branch  of  surface  and  the  rest  of  it  on  another  branch  of  surface. 

These  two  branches  of  surface  must  have  consistent  unit  normal  vectors  jn. 

If  one  branch  with  its  vector  n  rigidly  attached  to  it  was  rotated  about 

the  junction  line  until  it  coincided  with  the  other  branch,  the  vectors  n 

of  the  two  branches  should  point  in  opposite  directions.  The  expansion 

function  ua  is  unity  at  the  junction  node  and  exists  on  each  triangle 
junc 

that  is  simultaneously  on  one  of  the  two  branches  of  surface  in  the  domain 

of  u,  and  attached  to  the  junction  node.  On  any  one  of  these  triangles, 
junc 

u.  decreases  linearly  from  unity  at  the  vertex  which  is  the  junction  node 
junc 

to  zero  at  the  side  opposite  this  vertex.  The  expansion  function  for 

associated  with  u,  is  called  J.  .  According  to  (80),  J.  is  given 
junc  -junc  —junc 


J.  (r)  =  n  *  V  u,  (r) 
—junc  —  —  s  junc  — 


(105) 


As  many  such  linearly  independent  expansion  functions  as  possible  are 
associated  with  each  junction  node. 

Previously,  ST  separated  into  several  isolated  surfaces  when  some 
electrical  contacts  were  ignored  along  function  lines.  In  effect,  the 


L 


37 


function  u  in  (80)  was  specified  to  be  zero  at  several  points ,  one  on 

each  Isolated  surface.  When  the  electrical  contacts  are  restored, 

two  of  these  isolated  surfaces  may  combine  to  form  a  composite  surface. 

As  a  result, the  u's  on  both  of  these  formerly  isolated  surfaces  will 

combine  to  form  a  composite  u.  The  composite  u  is  not  necessarily 

continuous  on  the  junction  line  along  which  the  two  formerly  isolated 

surfaces  are  connected  to  each  other.  Presumably,  any  discontinuity  in 

u  gives  rise  to  an  impulse  in  V  u.  This  impulse  is  not  wanted  because, 

s 

according  to  (80),  it  would  appear  in  the  electric  current.  Fortunately, 
both  the  discontinuity  in  u  and  the  accompanying  impulse  in  Vgu  will  be 
suppressed  automatically  if  u  is  viewed  as  a  linear  combination  of  ex¬ 
pansion  functions,  each  of  which  exists  in  the  absence  of  all  the  others. 

We  now  generalize  (98)  to  Include  all  necessary  expansion  func¬ 
tions,  the  {uj}  of  (99),  the  in  (103),  the  in  (104),  and 

the  tujunc^  *n  (105).  It  is  more  convenient  to  test  (83)  before  inserting 
the  generalized  expansion  (98)  for  the  unknown  u.  Following  Galerkln's 
method  in  which  the  set  of  testing  functions  is  the  same  as  the  set  of 
expansion  functions,  we  choose  the  testing  functions  to  be  the  collection 
of  expansion  functions  in  the  generalized  expansion  (98).  This  collection 
of  functions  is  called  (u^,  i=l ,2, . . . ,Nl . 

The  integral  over  of  the  product  of  (83)  with  the  testing 
function  u^  is 


* 

J  . 


ui- 


V  x  A(J  )ds 


V  x  A(iabc)ds 


(106) 


We  want  u^ 


to  be  continuous  in  (106)  so  that  Stokes' 


theorem  can  be 


applied  to  both  integrals  in  (106).  Now,  may  be  of  (99),  may  be 
uhole  in  (103),  may  be  uloop  in  (104),  or  may  be  Ujunc  in  (105).  If 
is  Uj  of  (99)  or  u^0^e  in  (103),  then  is  continuous  on  S^,  the  range 
of  integration  in  (106). 

If  u.  is  u,  in  (104),  then  u.  drops  suddenly  from  unity  on  c. 

i  loop  i  J  loop 

to  zero  immediately  to  one  side  of  C^oQp.  In  this  case,  we  restrict  the 
range  of  Integration  in  (106)  to  the  part  of  ST  on  C^oop  and  to  the  side 
of  where  u^  decreases  linearly  from  unity  to  zero.  This  means  that 

the  part  of  to  the  other  side  of  C^oop  where  drops  suddenly  to  zero 

is  suppressed  from  the  range  of  integration  in  (106).  Not  affecting  the 


values  of  the  integrals  in  (106),  this  suppression  of  the  part  of  to 
one  side  of  c^oop  makes  Cloop  a  boundary  of  the  range  of  integration  so 
that  Uj  is  continuous  on  the  range  of  integration  and  equal  to  unity  on 
the  boundary  cloop* 

If  u,  is  u,  in  (105),  then  u.  exists  on  two  branches  of  surface 
i  June  i 

that  meet  on  a  junction  line.  Attaining  unity  on  the  junction  line,  u^ 
is  continuous  on  these  two  branches.  However,  an  observer  who  approaches 


the  junction  line  from  another  branch  of  surface  will  see  u^  jump  suddenly 
from  zero  on  that  branch  to  unity  on  the  junction  line.  Hence,  u^  is  dis¬ 
continuous  on  the  junction  line.  In  this  case,  we  restrict  the  range  of 
integration  in  (106)  to  the  two  branches  of  surface  on  which  u^^  exists. 

Not  affecting  the  values  of  the  integrals  in  (106),  this  restriction  of 
the  range  of  integration  renders  u^  continuous  in  (106).  Of  course,  this 
u^  is  zero  on  the  boundary  of  the  range  of  integration. 

With  the  range  of  integration  restricted  as  described  in  the  two 


previous  paragraphs,  (106)  becomes 
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11°  i* 


•  V  x  A(J  )ds 


■11“  i 2 


’  *  *<4bc)d» 


(107) 


where  S'  is  the  restricted  rsnge  of  integration.  Stokes'  theorem  con¬ 
verts  (107)  to  [3,  Eq.  (C-l)] 

|| (n  x  ysUi)  -AU^ds  -  J  u^C^)  *dj--jj  (nxV^)  ‘A^^ds 
c  •  r '  S' 


+ j  “ 


(108) 


where  C*  is  the  contour  that  bounds  S.  The  direction  of  C'  is  such  that 
a  right-handed  screw  would  advance  in  the  direction  of  n  when  rotated  in 
the  direction  of  C'.  Equation  (108)  is  rewritten  as 

|| ( nxV^)  ‘ACJ^ds  -  -  ||(]lxV8u1)  *A(Iabc)ds  +  |  t^ACJ)  *  dA  (109) 

S'  S'  c' 

where  J_  is  given  by  (20) . 

If  u^  is  Uj  of  (99)  or  Ujunc  in  (105),  then  u^  is  zero  on  C'  so 
that  (109)  reduces  to 

|  (nxV^)  •  AU^ds  -  -  ||  (n  x  V^)  •  A(Jabc)ds  (110) 

S'  S' 

If  u^  is  in  (103),  then  the  only  part  of  C'  on  which  u^  is  not  zero 

is  Chole.  On  Chole,  ut  is  unity  so  that  (109)  reduces  to 

|  (B  *  7sUi^  *A(Je)ds“-  |  (n  x  V^)  *A(Jabc)ds  +  |  A(J)  •  d l  (111) 

S’  8  S'  chole 

If  u  is  u,  in  (104),  then  the  only  part  of  C'  on  which  u.  is  not  zero 
i  loop  l 

is  C,  .  On  C,  ,  u.  is  unity  so  that  (109)  reduces  to 
loop  loop  i 

||  (n  x  V-Ul)  -A^ds  -  -  ||  (n  x  VuJ  -AU^ds  +  || 

cl  s'  c. 


A(J)  •  d£ 


(112) 


Both  and  C^oop  are  closed  contours.  In  view  of  (3),  Stokes' 

theorem  [2,  Eq.  (51)  on  p.  503]  gives 


and 


j  A(J)  •  d_£ 

cv  . 

hole 


!  B  (J) 
J 

^hole 


n  ds 


(113) 


J  A<i) 


c 


loop 


d  2 


j  B(J)  •  n  ds 
^loop 


(114) 


where  S,  ,  is  a  surface  that  caps  C.  ,  ,  and  S,  is  a  surface  that 
hole  r  hole  loop 

caps  c^oop-  Because  of  [3,  Eq.  (43)],  the  right-hand  sides  of  both 
(113)  and  (114)  vanish.  As  a  result,  both  (111)  and  (112)  reduce  to 
(110).  Therefore,  (109)  always  reduces  to  (110),  regardless  of  whether 
is  Uj  of  (99),  uhole  in  (103),  uloop  in  (104),  or  “junc  In  (105). 

If  all  the  impulses  in  V  u  are  suppressed,  (110)  can  be  written 

S  1 

as 


(n  *  V  u  )  •  A(J  )ds  =  - 
—  si - e 


(n  *  V  u 
s 


i>  •  iy,bc>ds 


(115) 


On  the  right-hand  side  of  (115),  the  integration  over  the  parts  of  S^,  on 

strips  a  and  c  is  difficult  to  perform  because  strips  a  and  c  on  which 

A(J  ,  )  must  be  evaluated  are  covered  with  the  source  current  J  .  .  Further- 

- abc  — abc 

more,  as  the  domain  of  u^  moves  on  either  strip  father  and  farther  from  the 
via,  the  right-hand  side  of  (115)  approaches  zero  in  a  manner  that  is  not 
obvious. 

Now,  the  objective  is  to  replace  the  right-hand  side  of  (115) 
by  a  form  that  is  more  suitable  for  calculation.  The  right-hand  side  of 


(115)  is  called  R  and  is  recast  as 
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R  -  R  +  R,  +  R 

a  1>  c 


where 


Ra  "  "  jj  *  W  *  A(4bc)  da 
S 

a 

*b  —  //  (*  x  W  *  ^(4bc)  ds 

sb 

Rc  "  -  j|  X  V8Ui)  *  -(4bC)  da 


(U6) 


(117) 


(118) 


(119) 


In  (117),  S  is  the  part  of  S  for  which  (-£.  £  z  <  z  ).  In  (118), 

A  X  A  A 

S.  is  the  part  of  S  for  which  (z  £  z  £  z  ) .  In  (119) ,  S  la  the  part 

of  S  for  which  (z  <  z  <  l  ) .  Otherwise  stated,  S  is  the  part  of  S 
Tec  A  1 

on  strip  a,  is  the  part  of  ST  on  the  via,  and  Sc  is  the  part  of  on 

strip  c.  As  in  (115),  all  impulses  are  suppressed  from  Vflu^  in  (117)-(119). 

Consider  Ra  of  (117).  According  to  (21),  £flbc  is  equal  to 

on  S  .  The  discontinuities  in  J  due  to  the  expansion  (84)  could  seriously 
A  —’A 

affect  the  accuracy  of  A(Jabc)  on  Sa»  Solving  (23)  for  and  substitu¬ 


ting  this  J  into  (26) ,  we  obtain 

— A 


j  ,  =  j.  +  j+  -  j+  +  j 

— abc  — b  — c  — a  — a 


Because  the  operator  A  is  linear,  (120)  allows  us  to  write 


A(J  .  )  -  A(J,  +  J+  -  J+)  +  A(J  ) 
—  — abc  -  — b  — c  —a a 


(120) 


(121) 


Substitution  of  (121)  into  (117)  gives 
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Ra  *  -  Jj  <H  *  Vi5  *  A<4+4"4)ds  '  j|  (-  X  Vsui)  *  (^(4)  ds  (122) 

So  S 

a  a 

Concerning  the  second  integral  in  (122),  Stokes'  theorem  yields 

[3,  Eq.  (C-l)] 

|  (n  x  V  u  )  •  A(J  )ds  ■  -  [  u  n  •  V x a(J  )ds  +  [  u  A(J  )  •  d£  (123) 

)  J  ~  si - a  J  J  i—  - a  J  i a  — 


where  C  is  the  contour  that  bounds  S  .  For  convenience,  S  is  extended 
a  a  a 

to  include  the  line  at  (z=z  ).  Now,  C  consists  of  the  lateral  edges  of 

a  a 

strip  a  for  (-i.  <  z  <  z  ),  the  edge  at  (z  ■  -  £,  ),  and  the  edge  at 

a  —  —  a  a 

(z  ”  z  ).  Because  u.  vanishes  on  the  lateral  edges  of  strip  a  and  because 
a  i 

A(J  )  is  z  directed  on  the  edge  at  (z  *  -<t  )  and  the  edge  at  (z  »  z  ),  the 
- a  a  a 

line  integral  over  C  on  the  right-hand  side  of  (123)  vanishes.  Moreover, 

& 

thanks  to  (40)  and  (41),  the  surface  integral  over  Sg  on  the  right-hand 
side  of  (123)  vanishes.  Hence,  (123)  reduces  to 


(n  x  y  u  )  •  A(J  )  ds  *  0 
—  si - a 


(124) 


so  that  (122)  becomes 


R  -  -  [  (n  x  V  u.)  •  A(J.  +  J+  -  J+)  ds  (125) 

a  J  j  —  si - b  — c  —a 


Consider  Rc  of  (119).  According  to  (21),  4bc  *s  e<lua^  t0  4 

on  S  .  The  discontinuities  in  J  due  to  the  expansion  (92)  could 
c  —c 

seriously  affect  the  accuracy  of  A(Jfl^c)  on  Solving  (24)  for 

and  substituting  this  J+  into  (26),  we  obtain 
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J  .  -  J.  +  J  -  J  +  J 
**abc  -i)  —a  — c  —c 


Because  the  operator  A  is  linear,  (126)  allows  us  to  write 


(126) 


A<4bc>  “  +  £  -  £>  +  A(4) 


(127) 


Substitution  of  (127)  into  (llg)  gives 


\  “  -  JJ  <H  x  Vi}  *  A(4+V4)ds  -  jj  (n  x  V^)  •  A(Jc)ds  (128) 

S  S 

c  c 

Concerning  the  second  integral  in  (128),  Stokes'  theorem  yields 

[3.  Eq.  (C-l)] 

JJ  (“  x  4V  *  A(Jc)ds  -  -  JJ  Uln  •  V  x  A(Jc)ds  +  J  u^AC^)  •  d£  (129) 


where  C£  is  the  contour  that  bounds  Sc>  The  right-hand  side  of  (129) 
vanishes  in  the  same  way  that  the  right-hand  side  of  (123)  did  so  that 


(128)  reduces  to 


R  *  -  [  (n  x  V  u.)  •  A(J,  +  J  -  J  )ds 

c  J  j  —  si  —  -4)  —a  — c 


(130) 


Substituting  (125),  (118),  and  (130)  into  (116),  we  obtain 

R  -  -  JJ  (n  x  VsUi)  •  A(4+4-4)ds  -  JJ  (n  x  7  uj  •  AU^ds 


[  (n  x  V  u.)  •  A(J.+J  -J  )ds 
JJ  —  si  —  — b  —a  — c 


(131) 


Hence,  (115)  becomes 
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j{  (£x  Vl> 

ST 

*  A(Je)ds  =  -  ||  (n 
ST 

x  V  u  )  •  A(J  ,  )ds 
s  i  —  — abc 

(132) 

Where  A(J  )  is  given  by 

(82),  and  A(Jabc)  is 

given  by 

f  Aub  + 

on  S 

a 

^(-abc)  = 

< 

on  S, 

D 

(133) 

/  A(4  +  Ja"  -  4) 

on  S 

c 

From  (101),  (103),  (104),  and  (105),  we  obtain 


J.  =  n  x  V  u. 
— i  -  si 


(134) 


where  all  the  impulses  in  V  u.  are  suppressed.  Equation  (134)  reduces 

S  1 

(132)  to 


A(J  )ds 


A  (J  .  ) ds 

- a  be 


(135) 


In  (135),  the  part  of  A(J  ^)  of  (133)  on  Sa  is  easy  to  evaluate 
because  the  points  of  evaluation  are  free  of  source  current.  Similarly, 
the  part  of  A (J^^)  °f  (133)  on  Sc  is  easy  to  evaluate  because  the  points 
of  evaluation  are  free  of  source  current.  Furthermore,  the  right-hand 


side  of  (135)  obviously  approaches  zero  as  the  domain  of  moves  farther 
and  farther  from  the  via  because  A  is  being  evaluated  farther  and  farther 
from  any  source  current.  On  S^,  however,  A  must  be  evaluated  on  the  source 
current  J[a^c •  According  to  (21),  J.aj5C  is  equal  to  on  S^.  We  choose 
to  be  a  constant  vector  on  each  triangular  patch  and  approximate  the  inte¬ 
gration  over  on  the  right-hand  side  of  (135)  by  sampling  A  at  the 
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centroids  of  the  triangles.  Since  the  centroids  of  the  triangles  are 
relatively  far  from  the  discontinuities  of  J_^  on  the  sides  of  the  tri¬ 
angles,  good  accuracy  should  be  obtained. 

In  view  of  (101),  (103),  (104),  and  (105),  the  generalized  ex¬ 
pansion  (98)  for  u  is  accompanied  by 


^  1=1  J  J 


(136) 


If  j  is  such  that  u^  is  one  of  the  uj's  1°  (98)  proper,  then  J_  is  given 

by  (102).  If  j  is  such  that  u^  is  u^o^e  1°  (103),  u^oop  in  (104),  or  u^  nc 

in  (105),  then  (102)  still  holds  provided  that  {T^,  n=l,2 . N_. }  are 

the  triangles  on  which  u.  exists,  A.  is  the  area  of  T,  ,  and  l,  is  an 

J  jn  jn  -jn 

appropriate  vector  side  of  T^.  As  defined  by  this  generalization  of  (102), 

the  generalized  flows  continuously  across  all  sides  of  triangles  in  the 

sense  that,  on  any  one  of  these  sides,  the  component  of  current  normal  to 

this  side  is  continuous.  Obviously,  the  generalized  J_  is  tangent  to  all 

triangle  sides  that  are  edges  of  ST> 

If  one  side  of  a  triangle  is  an  edge  of  ST>  then,  on  this  triangle, 

J  of  (136)  is  a  constant  vector  parallel  to  this  side.  However,  if  none 
— e 

of  the  sides  of  a  triangle  is  an  edge  of  S^,  then,  on  this  triangle, 
of  (136)  is  an  arbitrarily  directed  constant  vector.  Therefore,  of 
(136)  can,  if  necessary,  annihilate  on  any  triangle  on  the  via. 

Substituting  the  generalized  expansion  (136)  for  in  (135),  di¬ 
viding  (135)  by  I,  and  letting  (i=l , 2 , . . . ,N) ,  we  obtain  the  matrix  equation 


PI  -  V 


(137) 
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where  I  and  V  are  N  x 1  column  vectors,  and  P  is  a  square  matrix  of 
order  N.  The  jth  element  of  I  is  Ij  of  (136).  The  ith  element  of 
V  is  called  and  is  given  by 


J.  •  A(J  ,  )ds  ,  i=l,2,...,N 

—l - abc 


(138) 


The  ijth  element  of  P  is  called  P  and  is  given  by 


\ 


i=l ,2, . . .  ,N 


Pii  = 


J  •  A(J . )ds  ,  J 
— i - .1 


(j-1  ,2,  .  ...N 


Replacing  jn  by  im  in  the  generalization  of  (102),  we  obtain 


(139) 


J.  = 


5L 

— xm 


i  2A 


im 


on  T  ,  m=l ,  2 . N 

im  l 


(140) 


Substitution  of  (140)  for  J.  in  (138)  and  (139)  yields 


— l 


Ni  *. 

V  =  -  y 

i  L,  2IA.  J 

m=l  im  i 


A(J  ,  )ds ,  i=l ,2, . . . ,N 
—  —abc 


im 


(141) 


and 


P.  .  =  l 


i  £ .  f  r 

— im  !  ! 


i=l ,2,...,N 


ij  **-  2A. 

m=l  xm 


A(J.)ds  , 

T.  /  j=l  , 2  , .  .  .  ,N 

im 


(142) 


If,  in  the  integrands  of  (141)  and  (142),  A(^Ja^c)  and  A(J^) 
are  approximated  by  their  values  at  the  centroid  of  T^,  then  (141) 
and  (142)  reduce  to 


\--Ti  i ,4.  •  'A(i,bc>!Cl"  •  ‘■'•2 . N 

m*l 


(U3) 
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N. 
1  r1 


5  I  [A(J.)lclm  , 

2m=l~lm  --J  J-1.2 . N 


i=l , 2 ,  .  . .  ,N 


(144) 


where  the  superscript  cim  denotes  evaluation  at  the  centroid  of  T.  . 

Im 

Consider  the  A  that  appears  in  (144).  Since  A  is  given  by  (82) 
and  is  given  by  the  generalization  of  (102), 


[A(J  )]‘ 


N  . 

y  —1—  9  - - +  (u  t  -u  (t  -U  i  )  - — - 

?TT  n=i  Ajn  L_jn  Jj  |rCim-r*|  ~X  jnX  ~z  jnz  JJ  |rCim-r"|  - 


T 

jn 


(145) 


where  (£.  ,2.  ,2.  )  are  the  (x,y,z)  components  of  £.  ,  r'  is  the 

jnx  jny  jnz  v  '  — jn  — 

position  vector  of  the  differential  element  of  area  ds ' ,  r"  is  the  image 

of  jr'  about  the  ground  plane,  and  rclm  is  the  radius  vector  to  the  centroid 

of  T .  .  Each  integral  in  (145)  is  the  integral  over  a  triangle  of  the  re- 
lm 

ciprocal  of  the  distance  from  a  fixed  point.  The  value  of  this  integral 

is  given  by  [5,  Eq.  (46)],  an  expression  which  was  adapted  from  [9,  Eq .  (5)]. 

If  T.  is  on  either  S  or  S  ,  then  the  term  in  square  brackets  on  the  right- 
jn  a  c 

hand  side  of  (145)  has  no  y  component.  If  is  on  either  or  ,  then 

2  in  (144)  has  no  y  component  so  that  the  y  component  of  (145)  does  not 
— im 

come  into  play. 

Consider  [A(J  ,  )]Cim/I  of  (143).  Substituting  (26)  into  (133)  and 
—  —a  be 

using  the  fact  that  the  operator  A  is  linear,  we  obtain 

flACy]01"1/1  +  (A(J^)]Clm/I  -  lA(J^)lCinVl,  -^*<*3 

[A(Jabc)]cim/I  =<  [A(Jb)]clm/I  +  [A(j“)]Clm/I  +  [A(/)]Clm/I,  z^-z±zc 
/  [A(J  )]cln7l  +  [A(jj]Cim/I  -  [A(jj]clm/I,  z  <Z-?c 


(146) 
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where  the  superscript  cim  denotes  evaluation  at  the  centroid  of  T,  . 

im 


Now,  J  is  expressed  as 
— b 


-b  "  2Ahn 


°n  Tbn  ’  n=1’2'--->NN+l 


(147) 


where  (T^,  n=l,2, . . . ,N^+^)  are  the  triangles  on  which  J.  exists,  A^  is 
the  area  of  T^,  and  is  a  vector  that  is  independent  of  position  on 

Since  was  chosen  to  be  proportional  to  I  in  the  paragraph  pre¬ 
ceding  (98),  i,  can  not  depend  on  I.  Hence,  h  is  t ruly a  constant,  constant 
—bn  -on 

with  respect  to  spatial  coordinates  and  with  respect  to  I.  Moreover, 

{£,  }  must  be  such  that  J,  exists  only  on  the  via,  is  solenoidal,  is 

— bn  — b 

tangent  to  the  edges  of  S^,,  and  provides  for  a  continuous  flow  of  current 
from  J  on  line  a  through  the  via  to  J  on  line  c.  Lines  a  and  c  are  de- 

d  C 

fined  in  the  paragraph  that  contains  (21).  Construction  of  is  described 
in  Appendix  B.  Replacing  in  (82)  by  Jb  of  (147)  and  replacing  _r  in  (82) 
by  £C^m,  we  obtain 


fA(Jb)]Clm/I 


N+l  . 
_y_  V  l 

4tt  L  2A 

n=l  bn 


- ^ - +  (u  l,  -u  i,  -u  K  )  ds' 

i  cim  — y  bny  — x  bnx  — z  bnz  j  — — - 

ll  “I  I  T,  |rcim-r" 

bn  - 


(148) 

where  (1,  ,  ,  L  )  are  the  (x,y,z)  components  of  .  The  integrals 

bnx  bny  bnz  —bn 

in  (148)  are  evaluated  in  the  same  manner  as  those  in  (145). 

Next,  we  consider  (A(J  )]C^m/I  in  (146).  Since  J  is  given  by  (25a), 
—  — a  3 

[A(J_)]C*m  has  only  a  z  component  which  is  given  by  the  right-hand  side  of 


(47)  with  the  upper  limit  on  z'  replaced  by  z  .  Hence, 

cL 


49 


x  +w 
a  a 


[A(j")fim/I 


U  y^-rr  f  dx  '  .1  (x  '  )  dz'  ( 
-z  4ttI  J  az 


/,  cim  ,.22 
i/(z  -z')  +o 


f,  cim  i  ■.  2  2 

(z  -z  )  +p 

2 


,,  .  cim  cim  cim.  ,,  .  ,  ,  cim 

Here,  (x  ,  y  ,  z  )  are  the  (x,y,z)  coordinates  of  £  .If 

the  integral  with  respect  to  z'  in  (149)  is  called  I  then  f 4 , 


Expression  (152)  is  recast  as 


which  reduces  to 


Substitution  of  (154)  into  (149)  gives 


(149) 


where 


/  ,  cim  ,.2  .  cim  .  .  2 

=  /  (x  -x')  +  (y  -h  ) 

1  a 


(150) 


/,  cim  ,.2  .  cim,,  .2 

V  (x  -x ')  +  (y  +h  ) 


(151) 


Formula  200.01] 


r  *  i 1  a 

i  +  V  ,  +  Pi 

I  „  *  1  im  In  ( - ) 

z2  x  j— X - j 

-  y  +  Y  y  +  .  r.  - 

'  '  '  2  cim 

,=  -a-z 


(152) 


I  .  -  1  im  In  (- 
z2  a-*” 


/  2  2 

-r'+vY  +  P9  Pi 

- 4)  +  2  In  (— ) 


_  ,-z  -z 


2  .  2 


-Y+  V  Y  +  P 


Y=-a-z 


(153) 


cim  ,  /7c  im  .2  2 

z  -z  +  V (z  -z  )  +  p„ 

I  =  In  ( - - - ) 

z  cim  ,  cim  .2  2 

z  -z  +  V (z  ~za)  +  P, 

ci  a  * 


(154) 


x  +w 
a  a 


fA(jj]Clm/l  =  u 
—  —a  ; 


h 

z  4ttI 


c  im 
z  -z 


f.  cim  .2  2 

+  v  ( z  -z  )  +  p9 


J  (x1)  In  (- 
az 


/.cim  .2  2 

-z  +  v  (z  -z  )  +  p. 

a  a  i 


(155) 
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If  (zCim-z  )  in  (155)  is  viewed  as  the  axial  distance  from  the 
a 

end  of  the  strip  of  current  J  ,  then  (155)  generalizes  to 

3 


X  +w 
a  a 


r.,T±N1cim/T  u 

—(—a  /T  — z  4ttT 


J  (x  ' )  In  (■ 
az 


cim  \  ,  /  -  cim  .2,  2 

-f(z  -z  )  +  (z  - z  )  +  p 

- 5 -  a  — )dx' 

cim  .  cim  .2  2 

+(z  -z  )  +  V  (z  -z  )  +  Pi 

o  3  1 


(156) 


The  current  J  that  appears  in  (156)  is  given  by  (61)  which  is 

3  Z 


recast  as 


J  (x ' )  =  —  L  I  J  (x  ' ) 
az  M  a  az 


(157) 


where  L^,  the  inductance  per  unit  length  of  strip  a,  is  given  by  (67) 
which  is 


a  x  +w 

a  a 


(158) 


J  (x')dx' 
az 


Since  J  (x ' )  is  given  by  (84),  (158)  becomes 
az 


4ttA  y  I  . 

a  A 


(159) 


where  A  is  given  by  (86),  and  I  .  is  the  jth  element  of  the  column 
a  a  j 

vector  I  that  satisfies  (88).  Substitution  of  (84)  into  (157)  yields 
a 


4-rrL  I  a 

J  (x')  =  - —  l  I  .P  (x'-x  .) 

az v  d  aj  a'  aj 


(160) 


Substituting  (160)  for  J  ,  (150)  for  p. ,  and  (151)  for  f-  in  (156)  and 

3  Z  L  z. 
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then  changing  the  variable  of  integration  in  (156)  to  eliminate  the  offset 

cim  .  _  , 
x  ,  we  obtain 


- a  — z  a 


Nfl  xa.j+rx 

l  Xaj  J 

j-i  ajv  [ci 


+(zClm-z  )+  /x2+(yclm4li  )2+(*cim-z  )2 
ln( - 2 - —  . g  .  -  -  ■«  ,)dx 


xaj-xCim  +(zcln-Za)+  /x2+(ycim-ha)2+(zcim-*a)2 


(161) 

The  integral  in  (161)  is  of  the  same  form  as  the  integral  in  [5,  Eq.  (C— 4) ] . 
The  value  of  the  integral  in  [5,  Eq.  (C-4) ]  is  given  by  [5,  Eqa.  (C-23) 
and  (C-24)] . 

If  the  inductance  per  unit  length  of  strip  c  is  called  L  ,  then, 

c 

in  analogy  with  (159), 


(162) 


47rAc  E  lci 
c  j-1  CJ 


where  A£  is  given  by  (94),  and  1^  is  the  jth  element  of  the  column  vector 
I  that  satisfies  (96).  Replacing  J1  and  its  associated  quantities  in 

C  3 

+ 

(161)  by  and  its  corresponding  associated  quantities,  we  obtain 


n  xc,j+rx 


[A(J*)]Cim/I-u  L 
—  — c  Z  I 


I 

J  v  — : 


7/  cim  v  ,  /  2,  .  cira,.  >2,  ,  cim  72 
+(z  -z  )+Vx+(y  +h  )  +(z  -z  ) 

In  ( - - - -  -  - -  - — )  dx 


1  cim  7/  cim  .../  2.,  cim  .  .2.,  cim  .2 

x  .-x  +(z  -z  )+yx  +(y  -h  )  +(z  -z  ) 

CJ  c  c  c 

(163) 


The  integral  in  (163)  is  evaluated  in  the  same  manner  as  that  in  (161). 

Equations  (148),  (161),  and  (163)  are  substituted  into  (146),  and 

then  (146)  is  substituted  into  (143)  to  obtain  V^.  If  T^n  is  on  either  S# 

or  S  ,  then  i.  in  (143)  has  no  y  component  so  that  the  y  component  of 
c  — im 

(148)  does  not  come  into  play.  The  resulting  V.  is,  as  expected,  independent 
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of  I.  Substitution  of  (145)  into  (144)  gives  P  .  After  V.  and  P  .  have 

ij  i  ij 

been  calculated,  the  moment  solution  of  (83)  for  u  is  completed  by  solving 
(137)  for  I  and  substituting  the  elements  of  I  into  expression  (98)  for  u. 


V.  CONCLUSION 


To  conclude,  we  calculate  the  excess  inductance  by  substituting 

the  magnetostatic  currents  of  Section  IV  in  (38)  and  by  numerically 

performing  the  surface  integration  explicit  in  (38).  Calculated  values 

of  the  excess  current  density  on  the  strips,  the  total  current  density 

(J  +  J,  )  on  the  via,  and  the  excess  inductance  L  are  given  for  a  specific 
— e  — b  e 


example . 


Replacing  S  by  ST  in  (38)  and  using  the  fact  that  the  operator  A  is 


linear,  we  obtain 


L  =1  +1  +  I,  +  z  L  -  zL 

e  a  c  b  aa  cc 


where  L  is  given  by  (66),  and  L  is  given  by 

3  C 


(164) 


L  =  -T 
c  I 


(165) 


The  I's  in  (164)  are  given  by 


I  =  \  [A(.J  ,  )  +  A(J  )]  •  J  ds 

a  ,2  —  -abc - e  —a 


ic  -4  Jj  (Ay,bc>  •  4 ds 

1  s 

C 


4  -“2  JJ  +  ^(4)]  •  4  ds 

4 


(166) 


(167) 


(168) 
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i 


where  S  ,  S,  ,  and  S  are  defined  immediately  after  (119).  Furthermore, 
a  b  c 

A(J  bc)  is  given  by  (133)  in  which  d.abc  is  given  by  (26). 

To  numerically  perform  the  integrations  in  (166)  and  (167),  we 
enumerate  the  triangles  on  S  by  (T  ,  i=l,2,...,N  }  and  those  on 

HI  L  cl 

Sc  by  (T^,  i=Ntb+l,  N^+2 ,  •  •  •  ,Ntc } .  The  intervening  triangles 

{T.,i=N  +1,  N  +2,...,N  ,}  are  on  S,  and  will  come  into  play  later  when 

l  t3  tB  tb  b 

J  is  calculated.  Using  (25a)  to  replace  J  by  .1  in  (166)  and  then 
— e  —a  —a 

/v 

approximating  (A(J  )+  A(J  ))  on  each  triangle  of  S  in  (166)  by  its 

abCl  "0  B 

value  at  the  centroid  of  the  triangle,  we  obtain 


-7  (i(4bc>  + 

I  1=1 


(169) 


where  the  superscript  ci  denotes  evaluation  at  the  centroid  of  T^. 

Using  (25d)  to  replace  by  in  (167)  and  then  approximating 

(A(J  ,  )  +  A(J  ))  on  each  triangle  of  S  by  its  value  at  the  centroid 

- abc  —  — e  c 

of  the  triangle,  we  obtain 


N 

ic  -  -j  f  i*<4bc>  +  *<4),ci  ’  ![  4ds 

C  I  i=N  ,+1  3DC  i1 


(170) 


Replacement  of  n  by  m  in  (147)  and  subsequent  substitution  of 
(147)  into  (168)  produce 


1  N+l  Au  ff 

I.  ~  7  l  77^  •  [A(J  )  +  A(J  )]ds  ( 

b  1  m=l  2Abm  JJ  -~ahc 
bm 

Approximating  the  integral  in  (171)  by  sampling  the  integrand  at  the 


(171) 


centroid  of  T, 


we  obtain 
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,  NN+1 

h-2i  bc)+^4>' 

m=l 


cbm 


(172) 


where  the  superscript  cbm  denotes  evaluation  at  the  centroid  of  T 


bm 


As  for  J  in  (169),  combining  (45)  and  (160)  with  x'  replaced 
—a 


by  x,  we  obtain 


4ttL  1  a 

J  (x)  =  u  - —  l  I  .P  (x-x  .) 

—a  -z  M  A  aj  a  aj 


(173) 


where  L  is  given  by  (159),  and  I  .  is  the  jth  element  of  the  column 
sl  3  j 

vector  I  that  satisfies  (88).  On  T.,  J  (x)  of  (173)  is  approximated 

Si  13 


by 


4ttL  I 

,  ,  ,  a  Tci 

J  (x)  =  u - I 

—a  -z  y  a 


(174) 


-ci , 


where  I  is  the  average  value  over  x  on  T.  of  the  sum  in  (173). 
a  l 


,ci 


x  N 

i-max  a 


(illOA  « 

-5 - — -  [  l  I  .P  (x-x  ) 

i  i  .  ,L1aja  a.T 

x  -x  .  J 1  j=l 

max  min  x 


]  dx 


(175) 


min 


In  (175),  x1.  is  the  minimum  value  of  x  on  T  ,  and  x  is  the  maximum 
mm  i  max 

value  of  x  on  T..  Substitution  of  (174)  for  J  (x)  in  (169)  gives 
i  -a 


4rrL  Nta  . 

I-  I.  *daCl  +  *.y.»cl 


a  yl 


(176) 


i=l 


where  A  is  the  area  of  T.,  and  A  and  A  are  the  z  components  of  A 
i  i  z  z 

and  A,  respectively. 

As  for  .1  in  (170),  we  have,  in  analogy  with  (173), 


J  (x)  =  u 


N 

4ttL  I  c 
c 


-z  y 


y  i  . p  (x-x  ,) 

cj  c  cj 


(177) 
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where  L  is  given  by  (162),  and  I  .  is  the  jth  element  of  the  column 
c  c  j 

vector  that  satisfies  (96).  On  T^,  J^x)  of  (177)  is  approximated 


AuL  I 

J  (x)  =  u  -  I 

— c  —  z  p  c 


(178) 


where  I  is  the  average  value  over  x  on  T\  of  the  sum  in  (177). 


x  N 

/•max  c 


£  =,  __ - - -  [V  j  .P  (x-x  )  ]  dx 

c  i  i  ,  c  i  cv  cl' 


x  -x  .  J  1  j=l 
max  mm  x  , 
min 


(179) 


In  (179),  x1.  is  the  minimum  value  of  x  on  T,,  and  x^  is  the  maximum 
min  i  max 

value  of  x  on  .  Substitution  of  (178)  for  J^(x)  in  (170)  gives 


4tI.  tc 


Vlf,  J  Vc  '\<4bc>  +  V4>' 


(180) 


i=Ntb+l 


/*  /*v 

where  A,  is  the  area  of  T.  ,  and  A  and  A  are  the  z  components  of  A 
i  i  z  z  r  — 

and  A,  respectively. 

Expression  (176)  is  recast  as 


4ttL  ta 

X  =  - £  l  A.Icl 

a  U  x*]  i  3 


I  ^  (iiz«abc)ici/i  +  (»2(4)]ci/D 


and  (180)  is  recast  as 


(181) 


4ttL  tc 

c  r  .  Tci 


I  =  l  V  UA  (Jahr}1  /X  +  lA  O  /l} 

c  y  ,  vi  .  i  i  c  z  — abc  z  — e 


^Ntb+1 


(182) 


Similarly,  (172)  is  recast  as 


-  i  1,4.'  ,,A(4bc),cl“’/I  +  i*(4>'cb‘»> 

m*l 


(183) 
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The  next  task  is  to  evaluate  the  A^  and  A  terms  in  (181 ) — (183 ) 

and  the  A  and  A  terms  in  (181)-(183).  The  A  term  in  (181)  and  (182) 
z  z 

is  given  by  the  z  component  of  (146)  with  the  superscript  cim  replaced 
by  ci.  The  A  term  in  (183)  is  given  by  (146)  with  the  superscript  cim 
replaced  by  cbm.  As  indicated  in  the  final  paragraph  of  Section  IV, 
the  terms  on  the  right-hand  side  of  (146)  are  given  by  (148),  (161)  , 
and  (163). 


The  A  term  in  (181)  and  (182)  and  the  A  term  in  (183)  are  due 
z  — 

to  the  current  J^.  These  terms  are  evaluated  by  first  obtaining  a 

suitable  expression  for  J  .  Since  J  is  given  by  (136)  in  which  J. 

e  e  j 

is  given  by  (102),  we  have 


J 

— e 


N 

l 

j-1 


N. 

J 


J  n*l 


A, 

Z2IL 


2A 


.1" 


(184) 


where  I ^  is  the  jth  element  of  the  column  vector  I  that  satisfies 

(137),  A.  is  the  area  of  T.  ,  and  9.,  is  an  appropriate  vector  side 
Jn  jn  -jn 


To  eliminate 

the  overlapping  that  may  occur  in  (184)  in  the  sense  that  several  dif- 


of  T,  .  In  (184),  A.  is  viewed  as  existing  only  on  T 
jn  — jn  jn 


ferent  T,  's,  different  in  that  the  jn's  are  different,  may  actually 
jn 

be  the  same  physical  triangle,  we  rearrange  the  terms  in  (184)  into  a 


sum  of  vectors  each  of  which  exists  on  one  of  the  triangles  { T A } 

enumerated  in  the  paragraph  containing  (169).  Thus,  (184)  becomes 

N„ 

tc 

4  =  1  l  4 

^  k=i  ^ 


(185) 


where  is  a  constant  vector  on  T^.  Elsewhere,  4  zero*  0n 


4*  4  is  given  by 


I 
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I. 

.1  -jn. 


4  ‘  <  J  T 

J.n  jn  T 


(186) 


where  the  subscript  (T.  =  T.  )  denotes  that  Y  is  over  all  combi- 

jn  k  '• 

J  >n 

nations  of  j  and  n  in  (186)  such  that  (T,  =  T.  ) .  If  all  three 

jn  k 

vertices  of  are  on  an  edge  of  ST«  there  may  be  no  combination 

of  j  and  n  such  that  (T .  =  T.  )  in  which  case  Y  is  zero. 

in  k  “ 

1  >n 

Substituting  of  (185)  into  (82),  attaching  the  superscript 
ci  to  _r  and  A,  taking  the  z  component  of  (82),  and  dividing  by  t,  we 


obtain 


tc 

[A  (J  )]C1/I  =  •£-  l  K  ([[  - ^ - ff  — ^ - )  (187) 

k=l  k2  U  |  rC1-r  ’  |  [J  |rC1-r»l 


T  ■— 
k 


where  j:  *  is  the  radius  vector  to  the  differential  element  of  surface  ds', 

r"  is  the  image  of  _r'  about  the  ground  plane,  and  £  is  the  z  component 

of  As  in  (181)  and  (182),  the  superscript  ci  in  (187)  denotes  evalu- 

c  i 

ation  at  the  centroid  of  so  that,  in  (187),  _r  is  the  radius  vector 
to  the  centroid  of  T^.  Since  the  operator  A  is  given  by  (82)  and  the 


current  J  is  given  by  (185),  it  is  evident  that 


[A(JJ]cbm/I=^  l  (£, 


- - +(u  L  -u  £.  -u  JL  ) 

J]  |  cbm  , |  — y  ky  — x  kx  — z  kz  J J 

T,  -  1  T, 

k  k 


[  ds' 

Ji |rcbm_r"| 


(188) 


where  r',  ds',  and  r"  have  the  same  meanings  as  in  (187),  and 


^kx’Scy’^kz^  are  t*,e  components  of  As  in  (183),  the  super¬ 

script  cbm  in  (188)  denotes  evaluation  at  the  centroid  of  T  .  Thus,  in 


(188),  r 


.  The  integrals 


is  the  radius  vector  to  the  centroid  of  T, 
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in  (187)  and  (188)  are  evaluated  in  the  same  manner  as  those  in  (145). 

The  excess  inductance  L  can  now  be  calculated  from  (164)  with 

e 

1  given  by  (181),  I  by  (182),  I,  by  (183),  L  by  (159),  and  L  by  (162). 
a  c  b  a  c 

The  evaluation  of  the  A  and  A  terms  in  (181)-(183)  is  described  in  the 

z  — 

paragraph  that  follows  (183).  The  A^  term  in  (181)  and  (182)  is  given 

by  (187).  The  A  term  in  (183)  is  given  by  (188). 

A  computer  program  was  written  to  calculate  the  excess  inductance 

L  .  Intermediate  output  consists  of  the  excess  current  density  J  on  the 
e  — e 

strips  and  the  total  current  density  (J^  +  J^)  on  the  via.  This  computer 
program  will  be  described  and  listed  in  a  forthcoming  report.  Sample 
input  and  output  data  are  provided  to  verify  that  the  program  is  running 
properly.  For  the  sample  input  data,  N  =  8,  N  =10,  and  the  triangu- 
lated  surface  of  Fig.  4  is  used.  In  Fig.  4,  the  triangular  patches 

are  labeled  {T.,  j=l ,2 . 32  I .  The  vertices  of  the  triangles  are  called 

nodes  and  are  numbered  from  1  to  27.  All  dimensions  in  Fig.  4  are  in 
meters.  The  rise  from  (h^  =  3)  on  strip  a  to  (hc  =  5)  on  strip  c  is 
linear  on  the  via.  When  these  sample  input  data  were  entered,  the  com¬ 
puter  program  calculated  an  excess  inductance  of  0.704  micro-henrys 
and  the  current  density  shown  in  Tables  1  and  2.  The  permeability  U 
was  not  entered  as  data  but  was  indirectly  set  equal  to  the  permeability 
of  free  space  by  substituting  0.5  *  10~7  for  u  /  (8ti  )  near  the  end  of  the 


main  program. 
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Table  1.  The  x  component  of  current  density  J^/I  on  the  triangular 

patches  {T1+4(j_  i-1,2,3,4,  ,j=l  ,2, . . .  ,8  1  of  the  surface 

in  Fig.  4.  On  the  strips,  J  is  the  excess  current  density  J 

— e  — \ 

of  (185).  On  the  via,  J  is  the  total  current  density  J  +  J 

— e  3  —e  — 

where  J  is  given  by  (185)  and  J.  is  given  by  (B-l).  Table  1 

€t  D 

gives  the  x  component  of  J  /l  in  units  of  0.1  amperes  per 
meter . 


j  v 


1 

2 

3 

4 

5 

6 

7 

8 

4 

0.00 

0.00 

0.00 

2.29 

1 .40 

0.00 

0.00 

0.00 

1  3 

0.08 

0.18 

0.23 

2.11 

1 .04 

0.26 

0.14 

0.05 

1  2 

0.08 

0.18 

0.23 

1.13 

1.23 

0.26 

0.14 

0.05 

1 

0.0 

0.00 

0.00 

1.12 

1 .10 

0.00 

0.00 

0.00 

Table  2.  The  orthogonal  component  of  the  current  density  J  /i  on  the 

— e 

triangular  patches  {  ( j— 1 ) *  i=i»2>3,4,  j=l,2 . 8}  of  the 

surface  in  Fig.  4  where  is  defined  in  the  caption  for 

Table  1.  The  orthogonal  component  of  J^/l  is  the  component 

of  J  /I  perpendicular  to  the  x  direction.  This  component 

is  in  the  u  direction  on  the  strips  and  in  the  (2u  +  3u  ) 

— z  — y  — z 

direction  on  the  via.  Table  2  gives  the  orthogonal  component 
of  i  /I  in  units  of  0.1  amperes  per  meter. 


1 

2 

3 

4 

j 

5 

6 

7 

8 

4 

0.08 

0.26 

0.26 

1 .65 

1 .01 

-0.19 

-0.19 

-0.05 

3 

0.00 

0.08 

0.49 

2.99 

1  .65 

-0.42 

-0.05 

0.00 

2 

0.00 

-0.08 

-0.49 

1.98 

1.85 

0.42 

0.05 

0.00 

1 

-0.08 

-0.26 

-0.26 

— 

2.01 

_ 

1.98 

0.19 

0.19 

0.05 
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APPENDIX  A 


In  Appendix  A,  it  is  shown  that  If  A(r)  is  a  vector  function 
whose  component  tangent  to  a  surface  S  Is  differentiable  and  If 


n  •  V  x  A(r)  -  0  ,  r  on  S 


where  n  is  a  unit  vector  normal  to  S>  then  A.  (r)  can  be  written 

—  -■  f  nn 


(A-l) 


where 


A.  (r)  ■  V  ¥(r)  ,  r  on  S 
—tan  —  s  —  — 


f(r)  -  Y^)  +  J  A(r ')  •  dr’ 


(A-2) 


(A-3) 


In  (A-2),  the  subscript  tan  denotes  the  component  tangent  to  S,  and 

V  f(r)  is  the  surface  gradient  of  ¥(r).  In  (A-3),  r  is  the  position 
8  —  — 

vector  of  an  arbitrary  point  on  S,  and  C  is  any  contour  on  S  from 

r  to  r. 

— o  — 

The  following  reasoning  is  used  to  show  that  (A-2)  is  true. 
According  to  [2,  Eq.  166  on  p.  497] ,  we  have 

B  *  7  *  ^  “  h^  ^  (h2A2)  "  ^  ^iV1  (/ 

where  (v^,  v^)  are  orthogonal  curvilinear  coordinates  on  S,  and 
(h^,h2)  are  the  corresponding  metrical  coefficients.  Moreover,  A^ 
is  the  component  of  A  in  the  direction  of  Increasing  Vj ,  and  A2  is 
the  component  of  A  in  the  direction  of  Increasing  v2«  Equations 
(A-4)  and  (A-l)  imply  that 


(A-4) 


62 


4 


DVl  (h2A2)  "  Hv2  *hlAl* 


(A-5) 


In  view  of  (A-5),  the  differential  form 


hJA1dv]  +  h  A  dv 


(A-6) 


is  exact.  Therefore,  there  is  a  scalar  function  Y  such  that 


h1  A,  =  ~ 
1  1  .'v^ 

h2A2  ~  TvT 


(A-7) 

(A-8) 


From  (A-7)  and  (A-8),  we  obtain  the  desired  result  (A-2)  in  which 


„  ,u/  .  1  M  .  1  W 

V  1  (r)  =  7—  —  u.  +  —  t—  u 
s  —  h^  3v1  —1  h>2  »v2  -2 


(A-9) 


In  (A-9),  and  are,  respectively,  the  unit  vectors  in  the  direc¬ 
tions  of  increasing  v^  and  v2* 

If  £  changes  by  the  infinitesimal  amount  d_r,  then  4' (jr )  changes 
by  the  infinitesimal  amount  (V  V(r))  •  d£  so  that 


4' ( r )  =  Y(r  )  + 
—  — o 


(V'nr1))  •  dr’ 

s  -  — 


(A-10) 


where  r  and  C  are  the  same  as  In  (A—  3 ) .  In  (A-10),  V'Y(r')  is 
— o  r  s  — 

the  surface  gradient  of  4*  with  respect  to  the  coordinates  of  _r' . 
The  result  (A-3)  follows  from  (A-10)  and  (A-2). 


4 
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APPENDIX  B 


In  Appendix  B,  the  current  is  constructed  from  the  speci¬ 


fications  in  the  paragraph  containing  (147). 

In  agreement  with  (147),  we  expand  J,  as 

— b 


L> 

J.  -  I  l  1,  4J,  . 

*  *  1  — b  1 


(B-l  ) 


where  ( J,  }  are  vector  functions  each  of  which  is  constant  on  each 
— b 

triangle  on  the  via,  and  {l^j)  are  unknown  constants  to  be  determined. 

Since  J.  is  solenoidal  and  tangent  to  the  edges  of  S_,  it  is  natural 
— b  T 

to  choose  J  to  be  solenoidal  and  tangent  to  the  edges  of  S_.  In 
— bj  1 

order  for  to  be  solenoidal,  the  component  of  normal  to  any 

common  side  of  two  triangles  on  the  via  must  be  continuous  across  this 

side.  Here,  the  word  continuous  is  used  loosely  because,  if  bends 

sharply  at  the  side  of  a  triangle,  the  component  of  normal  to  this 

side  must  change  direction  suddenly  in  order  to  remain  on  S^,.  Guided 

by  these  conditions,  we  choose  some  of  the  functions  {J,  ,}  to  be  similar 

— bj 

to  J,  of  (104).  As  opposed  to  J  of  (102)  and  J,  of  (104),  the 
—loop  — J  —loop 

functions  (J  I  must  terminate  abruptly  on  line  a  and  on  line  c.  bine 
-bj 

a  is  the  line  at  (z=z  )  where  the  via  connects  with  strip  a.  Line  c 

a 

is  the  line  at  (z»z  )  where  the  via  connects  with  strip  c. 
c 

Now,  .1  must  provide  for  a  continuous  flow  of  current  from  3 
— b  a 

on  line  a  through  the  via  to  J  on  line  c.  Therefore,  on  line  a,  we 


must  have 
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J  ’  n  =  .1  (x) 

—  b  —  a  a  z 


(B-2) 


where  appears  in  (45),  and  n^  is  the  unit  vector  that  is  tangent 
to  on  the  via  side  of  line  a  and  perpendicular  to  line  a.  The 
vector  ri  points  from  line  a  toward  the  via.  On  line  c,  we  must  have 

3 


J  •  n  =  J  (x) 

— b  — c  c  z 


(B— 3) 


where  J  appears  in  (74),  and  n  is  the  unit  vector  that  is  tangent 
c  z  — c 

to  S^,  on  the  via  side  of  line  c  and  perpendicular  to  line  c.  The 

vector  n  points  from  the  via  toward  line  c.. 

— c 

If  has  N  nodes  on  line  a,  then  line  a  consists  of  (N„  -1) 

T  y.a  4a 

sides  of  triangles.  On  any  one  of  these  sides,  .J  (x)  varies  continu- 

3  Z 

ously  while  .J,  of  ( B— 1 )  remains  constant  so  that  (B-2)  can  not  be 
— b 

satisfied  everywhere  on  this  side.  Setting  J,  •  n_  equal  to  the  aver- 

D  3 

age  value  of  J  (x)  over  the  ith  triangle  side  on  line  a  and  letting  i 
az 

run  from  one  to  (N  -1),  we  obtain 

X  3 


.1  (x)dx  ,  i=l ,  2 ,  . .  .,No  -1 

az  *-a 


4*i  b.)  -b.i  -a  bat  v+  -  j 

^  “bai  bai  x 


(B-4) 


where  x,  is  the  x  coordinate  of  the  beginning  of  the  ith  triangle 
bai 

side  on  line  a,  and  x,  is  the  x  coordinate  of  the  end  of  the  ith 

ba  l 

triangle  side  on  line  a.  The  subscript  bai  on  the  left-hand  side  of 
(B-4)  denotes  evaluation  on  the  ith  triangle  side  on  line  a. 


If  S  has  N  nodes  on  line  c,  then  line  c  consists  of  (N .  -1) 

T  '■  c  *  c 

sides  of  triangles.  On  any  one  of  these  sides,  i^  ^(x)  varies  con¬ 
tinuously  while  of  (B-l)  remains  constant  so  that  (B-3)  can  not  be 
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satisfied  everywhere  on  this  side.  Setting  •  n^  equal  to  the 

average  value  of  J  (x)  over  the  kth  trianglj  side  on  line  c  and 
cz 

letting  k  run  from  one  to  (N^-l),  we  obtain 


n.  xi 

b  fb 

1  l  I.  .(J.  .  •  n  ),  ,  =  —r - — - 

. bj  — bi  — c  bck  +  -  I 

bck  bck  x. 


J  (x)dx,  k=l,2 . N„  -1 

c  z  ioc 


(B-5) 


where  x  is  the  x  coordinate  of  the  beginning  of  the  kth  tri- 
bck 

angle  side  on  line  c,  and  ,  is  the  x  coordinate  of  the  end  of 

bck 

the  kth  triangle  side  on  line  c.  The  subscript  bck  on  the  left-hand 
side  of  (B-5)  denotes  evaluation  on  the  kth  triangle  side  on  line  c* 

Together,  (B-4)  and  (B-5)  give  (N„  +  N.  -  2)  equations.  How- 

x.g  X,c 

ever,  one  of  these  equations  is  redundant  because,  due  to  the  nature 

of  {J,  the  flux  of  J,  entering  the  via  at  line  a  is  equal  to  the 

— bj  — b 

flux  of  J,  leaving  the  via  at  line  c,  and,  accordingly,  the  flux  of 
— b 

u  J  (x)  entering  the  via  at  line  a  is  equal  to  the  flux  of  u  J  (x) 
— z  az  — z  cz 

leaving  the  via  at  line  c.  Here,  u  is  the  unit  vector  in  the  z 
direction.  Now,  (N?a  +  -  3)  equations  remain  to  be  satisfied. 

Taking  the  number  of  functions  {J^ ^ }  in  ( B— 1 )  equal  to  the 
number  of  equations  to  be  satisfied,  we  have 


Nb  '  N«a  +  S«c  -  3 


<B-6> 


Each  of  the  first  (N„  -  1)  of  the  functions  {J.  . }  is  chosen  to  be 

y.a  - bj 

similar  to  J,  of  (104).  There  is  a  one-to-one  correspondence  be- 
—1  oop 

tween  these  (N^  -  1)  functions  and  the  (N^  -  1)  triangle  sides  on 
line  a.  The  function  associated  with  a  particular  triangle  side  on 
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line  a  is  not  zero  on  that  triangle  side  and  is  zero  on  all  other 

triangle  sides  on  line  a.  This  function  is  a  current  that  flows 

from  its  associated  triangle  side  on  line  a  through  the  via  to 

line  c..  Each  of  the  remaining  (N ,  -  2)  of  the  functions  {J.  . }  is 

v  e  — b.) 

chosen  to  be  similar  to  J  of  (102)  proper.  There  is  a  one-to-one 
correspondence  between  these  (N(,  -  2)  functions  and  the  (N^  -  2) 

interior  nodes  of  on  line  c.  The  function  associated  with  a 
particulai  interior  node  of  S  on  line  c  is  a  current  that  circulates 

L 

halfway  around  the  node.  Being  restricted  to  the  via,  this  current 
can  not  completely  encircle  the  node  as  I.  of  (102)  does. 

Dividing  both  sides  of  (B-4)  hy  I,  taking  all  (N  -  1)  equations 

*  a 

in  (B-4),  dividing  both  sides  ot  (B-5)  by  1,  and  taking  only  the  first 
(Nyc  -  2)  equations  in  (B-b) ,  we  obtain  the  matrix  equation 


Pi  -  V 
b  b  h 


(B— 7 ) 


where  V  and  T,  are  column  vectors,  and  P  is  a  square  matrix.  The  1th 
h  b  b 

element  of  I.  is  the  unknown  coeillclent  i,  .  In  (R-l).  The  ith  element 
b  b  ] 

of  V,  is  called  V,  .  and  is  given  by 
b  bi 


1 


bi  ,  +  w 

(x,  -x.  ,)T 

bat  bai 


bai 


ha  i 


J  (x)dx, 
az 


i“ 1  ,2, . . . ,N pfl-l 


(B-8) 


and 


(xbck”*bck)T 


hck 


bck 


I  (  X  )(lx  , 
i:  Z 


1  '  "«a'  "«*+1 . Nb  (B'5> 
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where  is  given  by  (B-6),  and  k  is  given  by 


k  =  i  -  N,  +  1 


(B-10) 


The  ijth  element  of  P^  is  called  P^^  and  is  given  by 


P.  .  .  =  (J.  . 
bij  -b.i 


and 


n  )  . 

—a  ba  1 


1-1,2 . N£a  1 


(B-ll ) 


1  Nia*  Nea+1,-,,,Nb 


P.  ..  = 

bij  -b] 


n  )u  I 

-c  bck 


(B-12) 


j=l  ,2  , .  . .  ,N, 


where  N  is  given  by  (B-6),  and  k  is  given  by  (B-10).  In  (B-ll),  n 
d  —a 

is  the  unit  vector  that  is  tangent  to  on  the  via  side  of  line  a  and 

perpendicular  to  line  a.  The  vector  n  points  from  line  a  toward  the 

—a 

via.  The  subscript  bai  in  (B-ll)  denotes  evaluation  on  the  ith  tri¬ 
angle  side  on  line  a.  In  (B-12),  n^  is  the  unit  vector  that  is  tangent 
to  on  the  via  side  of  line  c  and  perpendicular  to  line  c.  The  vector 
n^  points  from  the  via  toward  line  c.  The  subscript  bck  in  (B-12)  de¬ 
notes  evaluation  on  the  kth  triangle  side  on  line  c. 

Substitution  of  (61)  for  .T  (x)  in  (B-8)  gives 

az 


rbai 


bi 


x  +w 
a  r  a 


J  (x)dx,  i=l ,2, . . .  ,N0  -1 

n  Z  o 


(xbai_Xbai) 


.1  (x')dx'  bai 

az 


(B-13) 


where  J  (x)  satisfies  (57),  and  x,  .  and  x.  ,  are,  respectively,  the 
az  bai  bai 

x  coordinates  of  the  beginning  and  end  of  the  ith  triangle  side  on 


line  a.  Substitution  of  (75)  for  J  (x)  in  (B-9)  gives 

c  z 
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bi 


(xbck"Xbck) 


x  +w 
C  t  c 


bck 


JczMdx.  i-Nia,Nu+l,. 


■  tN 


J  (x')dx'  bck 
cz 


b 

(B-14) 


where  J  (x)  satisfies  (76),  N,  is  given  by  (B-6),  k  is  given  by  (B-10) , 
cz  b 

and  x,  ,  and  x*  ,  are,  respectively,  the  x  coordinates  of  the  beginning 
bck  bck 

and  end  of  the  kth  triangle  side  on  line  c. 

After  the  elements  of  P,  have  been  calculated  from  (B-ll)  and  (B-12) 

b 

and  those  of  Vt  from  (B-13)  and  (B-14),  the  matrix  equation  (B-7)  is 
b 

- ► 

solved  for  I,.  Then,  J,  is  given  by  (B-l)  in  which  I,  .  is  the  jth  ele- 
b  —  b  bj 


ment  of  I.  . 

b 
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